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Introduction 

Having  a  narrow  tibia  relative  to  body  mass  has  been  shown  to  be  a  major  predictor  of  stress  fracture 
risk  and  fragility  (Giladi  et  al,  1987;  Milgrom  et  al,  1989;  Beck  et  al,  1996).  The  reason  for  this 
phenomenon  is  not  understood.  Based  on  studies  of  genetically  distinct  inbred  mouse  strains,  we  found  a 
reciprocal  relationship  between  bone  mass  and  bone  quality,  such  that  slender  bones  are  associated  with 
more  damageable  bone  tissue  (Jepsen  et  al,  2001).  We  postulate  that  a  similar  reciprocal  relationship 
between  bone  mass  and  bone  material  properties  exists  in  the  human  skeleton.  The  intriguing  possibility 
that  slender  bones,  like  those  we  have  demonstrated  in  animal  models,  may  be  composed  of  more 
damageable  material  than  larger  bones  has  not  been  considered.  To  test  this  hypothesis,  we  propose  to 
determine  whether  whole  bone  geometry  is  a  predictor  of  tissue  fragility  in  the  tibiae  from  young  male 
donors.  Tissue  damageability  will  be  assessed  from  biomechanical  testing  of  compact  bone  samples  and 
correlated  with  measures  of  bone  slenderness.  Specimens  will  be  subjected  to  detailed  analyses  of  bone 
microstructure,  composition,  and  microdamage  content.  In  the  second  set  of  experiments,  these  analyses 
will  be  repeated  for  female  donors  to  test  for  gender  differences  in  tissue  fragility. 

Body 

In  the  fourth  year  of  this  grant,  we  focused  on  continuing  the  tissue-level  mechanical  tests  for  the 
female  bones  and  the  histological  analysis  of  the  male  bones.  We  applied  for  a  no-cost  extension  to 
continue  this  research  in  order  to  finish  the  last  half  dozen  female  samples.  Obtaining  human  samples 
from  females  ages  18-40  has  been  a  challenge  recently  and  we  have  expanded  our  resources  to  include 
the  Anatomy  Gifts  Registry  and  the  National  Disease  Research  Interchange. 

In  year  three,  we  had  begun  a  series  of  experiments  that  examined  whether  ultrasound  has  the 
sensitivity  to  detect  the  presence  of  damage  in  bone.  Most  studies  have  found  that  ultrasound  is 
insensitive  to  damage.  We  have  found  that  damage  affects  the  viscoelastic  properties  of  bone  to  a  greater 
degree  than  the  stiffness  properties  and  we  have  focused  our  attention  on  the  attenuation  of  the 
ultrasound  signal  as  this  may  be  more  sensitive  to  the  presence  of  damage.  A  student  for  the  City 
University  of  New  York  worked  on  this  project  and  found  that,  like  prior  work,  the  velocity  of  the 
ultrasound  signal  was  insensitive  to  the  presence  of  damage.  He  examined  the  attenuation  (a  measure 
related  to  the  viscoelastic  properties  of  bone)  following  damage  accumulation,  however,  could  not  make 
a  reasonable  conclusion  because  of  difficulties  in  matching  the  size  of  the  sample  to  the  size  of  the 
ultrasound  transducer.  Although  viscous  properties  may  be  more  sensitive  to  the  presence  of  damage,  we 
concluded  from  this  study  that  it  was  not  practical  to  reliably  measure  this  parameter  using  machined 
bone  samples  from  human  tibiae.  We  would  have  to  develop  protocols  using  a  larger  animal  (e.g.,  cow) 
to  further  test  this. 

A  major  outcome  of  this  grant  was  finding  that  tissue-level  mechanical  properties  vary  with  the  size 
of  the  tibia  (Tomassini  et  al,  JBMR,  2005).  Tissue  modulus  decreased  with  bone  size,  whereas  post¬ 
yield  strain  and  total  energy  increased  with  bone  size.  Not  all  of  these  correlations  were  significant 
indicating  that  the  relationship  between  bone  morphology  and  tissue  level  mechanical  properties  was 
subtle.  The  damage  parameter  decreased  with  bone  size  but  only  when  the  morphological  properties 
were  nonnalized  for  body  weight  and  size.  These  results  indicated  that  smaller  bones  are  comprised  of 
tissue  that  is  more  stiff  and  less  ductile  (i.e.,  more  brittle)  compared  to  larger  bones.  In  year  four,  we 
conducted  a  histological  evaluation  to  detennine  if  we  could  explain  the  variability  in  tissue-level 
mechanical  properties  based  on  differences  in  tissue  microstructure.  The  mid-span  region  of  the  samples 
that  were  subjected  to  damage  accumulation  tests  were  stained  with  basic  fuchsin  and  embedded, 
undecalcified,  in  poly-methylmethacrylate.  Transverse  sections  (150micron  thickness)  were  cut  using  a 
low-speed  diamond  waffering  saw,  ground  and  polished  to  lOOmicron  thickness,  and  mounted  on  glass 
slides  for  evaluation.  Digitial  images  of  each  cross-section  were  analyzed  for  the  area  fraction  of 
Haverisian  porosity,  the  area  fraction  of  non-Haverisian  porosity,  the  total  porosity,  the  area  fraction 
occupied  by  full  osteons,  the  area  fraction  occupied  by  osteon  fragments,  and  the  area  fraction  occupied 
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by  circumferential  (unremodeled)  bone.  Further,  an  adjacent  portion  of  the  test  samples  were  subjected 
to  standard  ashing  methods  in  order  to  assess  ash  content.  We  corrected  the  microstructural  and 
compositional  values  for  age-effects  by  linear  regression  analysis  and  then  constructed  a  correlation 
matrix.  We  found  that  the  post-yield  strain  (p<0.04)  and  total  energy  (p<0.06)  were  positively  correlated 
with  the  area  fraction  of  full  osteons.  No  other  significant  correlations  were  observed.  This  indicated  that 
having  a  lower  area  fraction  of  osteons  lead  to  a  more  brittle  failure  mode.  No  significant  correlations 
were  observed  between  ash  content  and  any  of  the  mechanical  properties  tested.  This  contrasts  with  the 
results  for  the  mouse  femurs  which  showed  that  variation  in  whole  bone  ductility  was  explained  largely 
by  ash  content  (Jepsen  et  al,  Mamm  Genome,  2004). 


Figure  1.  Correlation  between  a)  post-yield  strain  and  b)  total  energy  and  the  area  fraction  of  full 
osteons. 
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We  have  analyzed  12  of  the  18  female  tibia  to  date,  and  are  aggressively  trying  to  locate  6  more 
pairs  of  tibia.  The  problem  with  our  current  data  set  is  that  1 1  of  the  12  females  range  in  age  from  33-46. 
Consequently,  the  data  set  is  not  easily  related  to  the  males  until  we  obtain  additional  samples  in  the  20- 
30  year  old  range.  Obviously  this  is  a  challenge  since  this  age  group  and  gender  is  one  of  the  rarest  to 
find,  particularly  for  samples  that  are  free  of  bone-related  illnesses.  We  are  on  the  wait-list  of  several 
bone  transplant  organizations,  and  we  anticipate  that  we  should  be  able  to  secure  the  remaining  6 
samples  over  the  next  year. 
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Key  Research  Accomplishments 

The  primary  outcome  of  the  grant  thus  far  is  that  individuals  with  more  slender  bones  appear  to  have 
different  material  properties.  With  increasing  slenderness,  the  material  is  more  stiff,  less  ductile,  and 
more  damageable.  This  supports  are  central  hypothesis  and  may  help  explain  why  individuals  with 
smaller  tibiae  are  at  higher  risk  of  stress  fractures. 

Reportable  Outcomes 

Bouxien,  ML,  Jepsen  KJ.  Etiology  and  biomechanics  of  hip  and  vertebral  fractures.  Atlas  of 
Osteoporosis,  Second  Edition.  Current  Medicine,  Inc.,  Eds.  Eric  S.  Orwoll,  Stanley  G.  Korenman, 
2003. 

Jepsen  K.  The  aging  cortex:  to  crack  or  not  to  crack.  Osteoporos  hit.  2003  Sep;  14  Suppl  5:57-66.  2003. 
Tommasini  SM,  Nasser  P,  Jepsen  KJ.  Gender  differences  in  bone  slenderness  are  not  related  to  material 
properties.  Transactions  Orthopaedic  Research  Society,  2002. 

Tommasini  SM,  Morgan  TG,  van  der  Meulen  MCH,  Jepsen  KJ.  Genetic  variation  in  vertebral 
mechanical  properties  detennined  by  the  relationship  between  morphological  and  compositional 
bone  traits.  Transactions  Orthopaedic  Research  Society,  2003. 

Jepsen  KJ,  Price  C,  Nadeau  JH.  Systems  analysis  of  bone  fragility.  Pathways,  Networks,  and  Systems: 

Theory  and  Experiments.  Aegean  Conference,  2003. 

Tommasini  SM,  Nasser  P,  Jepsen  KJ.  The  relationship  between  bone  morphology  and  bone  quality: 
Implications  for  stress  fracture  risk  in  young  adult  male  tibiae.  Poster  and  podium  presentations  as 
the  American  Society  of  Bone  and  Mineral  Research  Annual  meeting,  Seattle,  WA,  2004. 

Bird  JE,  Nasser  P,  Tommasini  S,  Casagrande  D,  Jepsen  KJ.  The  relationship  between  continued 
periosteal  apposition  and  bone  fragility.  Poster  presentation  at  the  American  Society  of  Bone  and 
Mineral  Research  Annual  meeting,  Seattle,  WA,  2004. 

Tommasini  SM,  Nasser  P,  Schaffler  MB,  Jepsen  KJ.  The  relationship  between  bone  morphology  and 
bone  quality  in  male  tibiae:  Implications  for  stress  fracture  risk.  Journal  of  Bone  and  Mineral 
Research,  20(8):  1372-1380,  2005. 

Price  C,  Herman  BC,  Lufkin  T,  Goldman  HM,  Jepsen  KJ.  Genetic  variation  in  bone  growth  patterns 
defines  adult  mouse  bone  fragility.  Journal  of  Bone  and  Mineral  Research,  20(11):  1983-1991, 
2005. 

Jepsen,  KJ.  Assessing  the  mechanical  properties  of  mineralized  tissues.  Invited  talk,  Advances  in 
Mineral  Metabolism  (AIMM);  Snowmass,  CO,  2005. 

Jepsen,  KJ.  Using  mice  to  understand  structure-function  relationships  in  the  skeleton.  Invited  talk,  Bone 
Quality  Meeting;  Sponsored  by  the  National  Institutes  of  Arthritis  and  Musculoskeletal  and  Skin 
Diseases  (NIAMS)  and  the  American  Society  of  Bone  and  Mineral  Research;  Bethesda  MD  May  2- 
3,  2005. 

Jepsen,  KJ.  Biomechanical  insights  into  the  components  of  bone  strength.  Invited  talk,  Alliance  for 
Better  Bone  Health  Sponsored  Symposium  at  the  Second  Joint  Meeting  of  the  European  Calcified 
Tissue  Society  (ECTS)  and  the  International  Bone  and  Mineral  Society  (IBMS),  Geneva, 
Switzerland,  June,  2005. 

Jepsen,  KJ.  Can  heritable  aspects  of  bone  quality  be  enhanced  through  intervention?  Invited  talk, 
American  Society  of  Bone  and  Mineral  Research;  Working  group  on  Bone  Remodeling  and  Stress 
Fracture,  Nashville,  TN,  2005. 

Funding 

The  data  generated  by  this  grant  provided  evidence  that  the  mouse  represents  an  important  model  for 
understanding  the  genetic  variability  in  the  human  skeleton.  This  data  helped  secure  an  ROl  grant  from 
the  NIH  (AR44927)  titled,  "Genetic  Determination  of  Skeletal  Fragility". 
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The  data  generated  by  this  grant  was  also  used  in  the  submission  of  two  NIH  grants  titled,  “Genetic 
regulation  of  bone  growth  and  development”  and  “Age-changes  in  bone  morphology  as  a  predictor  of 
fracture  risk”.  The  goal  of  the  first  grant  is  to  identify  chromosomes  harboring  quantitative  trait  loci 
influencing  bone  growth  patterns.  The  goal  of  the  second  grant  is  to  detennine  if  25  year  changes  in 
moment  of  inertia  predict  fracture  incidence.  This  latter  grant  utilizes  information  from  the  Framingham 
Heart  Study  and  was  written  in  collaboration  with  the  Boston  University  School  of  Public  Health  and  the 
Hebrew  Rehabilitation  Center  for  the  Aged. 

Conclusions 

The  results  to  date  have  provided  new  insight  into  the  relationship  between  bone  morphology  and  tissue 
mechanical  properties.  The  investigations  of  the  mouse  skeleton  revealed  that  genetic  variation  in  bone 
morphology  strongly  influence  tissue  mechanical  properties  through  variations  in  matrix  composition. 
The  data  suggest  that  a  similar  relationship  may  also  exist  in  the  human  skeleton.  Thus,  individuals  who 
have  smaller  (more  narrow)  tibia  for  their  body  size  may  compensate  for  the  smaller  geometry  through 
variation  in  tissue-level  mechanical  properties.  One  of  the  side  effects  of  this  compensation  is  altered 
damageability  which  may  be  revealed  under  extreme  physical  activity  such  as  that  experience  during 
military  training. 
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Relationship  Between  Bone  Morphology  and  Bone  Quality  in  Male 
Tibias:  Implications  for  Stress  Fracture  Risk 


Steven  M  Tommasini,1  Philip  Nasser,2  Mitchell  B  Schaffler,2  and  Karl  J  Jepsen2 


ABSTRACT:  Biomechanical  properties  were  assessed  from  the  tibias  of  17  adult  males  17-46  years  of  age. 
Tissue-level  mechanical  properties  varied  with  bone  size.  Narrower  tibias  were  comprised  of  tissue  that  was 
more  brittle  and  more  prone  to  accumulating  damage  compared  with  tissue  from  wider  tibias. 

Introduction:  A  better  understanding  of  the  factors  contributing  to  stress  fractures  is  needed  to  identify  new 
prevention  strategies  that  will  reduce  fracture  incidence.  Having  a  narrow  (i.e.,  more  slender)  tibia  relative  to 
body  mass  has  been  shown  to  be  a  major  predictor  of  stress  fracture  risk  and  fragility  in  male  military  recruits 
and  male  athletes.  The  intriguing  possibility  that  slender  bones,  like  those  shown  in  animal  models,  may  be 
composed  of  more  damageable  material  has  not  been  considered  in  the  human  skeleton. 

Materials  and  Methods:  Polar  moment  of  inertia,  section  modulus,  and  antero-posterior  (AP)  and  medial- 
lateral  (ML)  widths  were  determined  for  tibial  diaphyses  from  17  male  donors  17-46  years  of  age.  A  slen¬ 
derness  index  was  defined  as  the  inverse  ratio  of  the  section  modulus  to  tibia  length  and  body  weight.  Eight 
prismatic  cortical  bone  samples  were  generated  from  each  tibia,  and  tissue-level  mechanical  properties  in¬ 
cluding  modulus,  strength,  total  energy,  postyield  strain,  and  tissue  damageability  were  measured  by  four- 
point  bending  from  monotonic  ( n  =  4/tibia)  and  damage  accumulation  (n  =  4/tibia)  test  methods.  Partial 
correlation  coefficients  were  determined  between  each  geometrical  parameter  and  each  tissue-level  mechani¬ 
cal  property  while  taking  age  into  consideration. 

Results:  Significant  correlations  were  observed  between  tibial  morphology  and  the  mechanical  properties  that 
characterized  tissue  brittleness  and  damageability.  Positive  correlations  were  observed  between  measures  of 
bone  size  (AP  width)  and  measures  of  tissue  ductility  (postyield  strain,  total  energy),  and  negative  correlations 
were  observed  between  bone  size  (moment  of  inertia,  section  modulus)  and  tissue  modulus. 

Conclusions:  The  correlation  analysis  suggested  that  bone  morphology  could  be  used  as  a  predictor  of  tissue 
fragility  and  stress  fracture  risk.  The  average  mechanical  properties  of  cortical  tissue  varied  as  a  function  of  the 
overall  size  of  the  bone.  Therefore,  under  extreme  loading  conditions  (e.g.,  military  training),  variation  in  bone 
quality  parameters  related  to  damageability  may  be  a  contributing  factor  to  the  increased  risk  of  stress  fracture 
for  individuals  with  more  slender  bones. 
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INTRODUCTION 

Stress  fractures  are  overuse  injuries  of  bone  that  are 
common  among  elite  runners  and  military  recruits/1-3) 
Before  injury,  affected  bones  are  typically  normal  with  no 
acute  injury.  Morbidity  from  stress  fractures  ranges  from 
minor  pain  to  serious  lifetime  disability  for  the  individual/4) 
Stress  fractures  have  been  reported  in  the  ribs,  hip,  spine, 
and  metatarsals/3’5^  but  vigorous  weight-bearing  activities, 
such  as  running  and  jogging,  commonly  lead  to  stress  frac¬ 
tures  of  the  lower  extremities,  especially  the  tibia/3)  During 
basic  training,  1-5%  of  U.S.  male  military  recruits  sustain  a 
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stress  fracture/2-1  However,  this  incidence  is  two  to  five 
times  higher  in  female  recruits/6)  Stress  fractures  lead  to 
loss  of  manpower,  valuable  loss  of  training  time,  expense  of 
medical  care,  and  discharge  of  affected  soldiers/71  A  better 
understanding  of  the  factors  contributing  to  stress  fractures 
is  needed  to  identify  new  prevention  strategies  that  will 
reduce  fracture  incidence. 

A  number  of  risk  factors  for  stress  fracture  have  been 
identified  including  physical  fitness,  external  hip  rotation, 
body  height  and  weight,  age,  race,  gender,  muscle  mass, 
motivation,  footwear,  smoking,  and  family  history  of  osteo¬ 
porosis/1, 4,8_1°)  One  of  the  best  predictors  of  stress  fracture 
risk  is  bone  geometry.  Specifically,  having  a  narrow  (i.e., 
more  slender)  tibia  relative  to  body  mass  has  been  shown  to 
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be  a  major  predictor  of  stress  fracture  risk  and  fragility  in 
male  military  recruits(1’2,11’  and  male  athletes/12’  A  stress 
fracture  is  thought  to  be  a  consequence  of  transiently  re¬ 
duced  tissue  strength  arising  from  increased  resorptive  ac¬ 
tivity  (i.e.,  increased  porosity)  that  acts  to  repair  damage 
induced  by  vigorous  physical  activity/13’  Thus,  stress  frac¬ 
tures  may  be  pronounced  in  individuals  with  more  slender 
bones  because  smaller  bone  size  is  thought  to  lead  to  higher 
tissue-level  stresses  and  thus  increased  damage  accumula¬ 
tion/1,2’  However,  this  postulate  is  based  on  the  assumption 
that  all  bones  are  constructed  in  equivalent  manners,  and 
the  contribution  of  variable  tissue-level  mechanical  proper¬ 
ties  to  stress  fracture  incidence  has  not  been  explored. 

An  examination  of  inbred  mouse  strains  may  help  ex¬ 
plain  why  bone  size  is  a  risk  factor  for  stress  fractures  in  the 
human  skeleton.  A  comparison  of  adult  A/J  and  C57BL/6J 
inbred  mouse  strains  revealed  that  the  bone  slenderness 
was  inversely  related  to  mineral  content  (as  measured  by 
ash  content)  and,  by  correlation,  tissue  modulus  and 
strength/14’  Mineral  content  has  been  shown  to  be  posi¬ 
tively  correlated  with  tissue  stiffness  and  strength/15’  These 
results  suggested  that  bone  morphology  and  mineral  con¬ 
tent  were  coordinately  regulated  so  whole  bone  stiffness 
appropriately  matched  the  mechanical  demands  imposed 
by  weight  bearing.  However,  the  downside  of  regulating 
mineral  content  to  match  bone  size  was  that  mineral  con¬ 
tent  was  also  negatively  correlated  with  tissue  ductil- 
ity/14,15’  We  postulate  that  a  similar  reciprocal  relationship 
between  bone  size  and  bone  quality  exists  in  the  human 
skeleton.  The  intriguing  possibility  that  slender  bones,  like 
those  shown  in  animal  models,  may  be  composed  of  more 
damageable  material  has  not  yet  been  considered  in  the 
human  skeleton. 

The  goal  of  this  study  was  to  determine  whether  tissue- 
level  mechanical  properties  vary  with  bone  size  in  the  human 
skeleton.  This  was  tested  by  assessing  the  biomechanical 
properties  of  tibias  from  young  adult  males.  Understand¬ 
ing  why  bone  morphology  is  a  risk  factor  for  stress  frac¬ 
tures  should  lead  to  better  identification  of  those  at  risk 
and,  ultimately,  to  early  diagnosis,  treatment,  and  modifi¬ 
cation  of  training  regimens. 

MATERIALS  AND  METHODS 

Sample  population 

Tibias  of  17  male  donors  (15  white,  1  Hispanic,  1  black) 
32.9  ±  10.4  years  of  age  (range,  17^16  years)  were  acquired 
from  the  Musculoskeletal  Transplant  Foundation  (Edison, 
NJ,  USA).  Donor  body  weight  and  height  were  obtained 
from  the  source.  Only  donors  with  no  known  skeletal  pa¬ 
thology  were  included  in  the  study.  The  tibias  were  freshly 
harvested,  wrapped  in  wet  gauze,  and  stored  in  plastic  bags 
at  -40°C. 

Whole  bone  morphology 

Tibia  length  (L)  was  measured  as  the  average  distance 
between  the  distal  articular  center  (the  middle  of  the  talar 
trochlear  facet)  and  the  two  proximal  articular  centers  (me¬ 
dial  and  lateral  condyles)(16’  using  a  large-capacity  slide 
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FIG.  1.  Schematic  of  how  whole  tibias  were  sectioned  to  pro¬ 
duce  the  3-mm-thick  sections  used  for  cross-sectional  morphology 
and  cortical  bone  samples  for  biomechanical  testing  from  three 
diaphyseal  cylindrical  sections  (monotonic,  n  =  4/tibia;  damage 
accumulation,  n  =  4/tibia). 


caliper  with  an  accuracy  of  ±2.54  mm  (Mantex  Precision; 
Haglof,  Madison,  MS,  USA).  Tibia  width  was  measured  in 
the  antero-posterior  (AP)  and  medial-lateral  (ML)  direc¬ 
tions  at  10%  intervals  from  30-70%  of  the  total  tibia  length 
using  a  300-mm  vernier  caliper  with  an  accuracy  of  ±0.02 
mm  (Fowler  Company,  Newton,  MA,  USA). 

Cross-sectional  morphology  was  determined  from  3-mm- 
thick  middiaphyseal  cross-sections  cut  at  30%,  50%,  and 
70%  of  the  total  tibia  length  (Fig.  1)  using  a  diamond 
coated  metallurgical  saw  (Model  660;  South  Bay  Technol¬ 
ogy,  San  Clemente,  CA,  USA).  A  calibrated  image  of  each 
cross-section  was  obtained  using  a  digital  camera  at  a  0.024 
mm/pixel  resolution.  Image  analysis  software  (IMAQ  Vi¬ 
sion  Builder  6.0;  National  Instruments,  Austin,  TX,  USA) 
was  used  to  threshold  each  image  and  quantify  cortical  area 
(CtAr),  the  moments  of  inertia  about  the  AP  (/AP)  and  ML 
(/ML)  axes,  the  polar  moment  of  inertia  (J  =  IAP  +  /ML), 
and  the  section  modulus  in  the  AP  (J/APwidth/2)  and  ML 
(J/MLwidth/2)  directions.  Moment  of  inertia  and  section 
modulus  were  assessed  because  these  geometric  measures 
are  related  to  the  bending  and  torsional  stiffness  of  intact 
tibias.  A  slenderness  index  (S)  was  calculated  in  the  AP  and 
ML  directions  as  the  ratio  of  the  AP  and  ML  section  modu¬ 
lus  values,  respectively,  to  tibia  length  and  body  weight(17’: 

S  =  l/{[J/(width/2)]/(L  x  BW))  (1) 

where  L  =  tibia  length  (mm)  and  BW  =  body  weight  (kg). 
The  section  modulus  has  been  shown  to  scale  linearly  with 
body  mass/17’  The  inverse  ratio  was  used  so  that  a  tibia 
with  a  large  slenderness  value  is  one  that  is  thinner  or  grac- 
ile  for  the  weight  and  height  of  an  individual.  A  small  slen¬ 
derness  value  reflects  a  stocky  or  robust  tibia.  All  morpho¬ 
logical  traits  were  averaged  over  the  three  cross-sections  for 
each  tibia. 

Bone  sample  generation 

Cortical  bone  samples  were  prepared  from  the  diaphysis 
of  each  tibia  for  biomechanical  testing  (Fig.  1).  The  three 
diaphyseal  cylindrical  sections  were  rough-cut  into  antero- 
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FIG.  2.  (A)  Schematic  of  the  four-point  bending  device  shows 

how  cortical  bone  samples  were  tested.  The  span  between  the 
upper  contact  points  is  2a  =  18  mm,  and  the  span  between  the 
lower  contact  points  is  2L  =  42  mm.  (B)  Typical  stress-strain 
curve  from  the  four-point  bending  monotonic  tests  shows  how 
modulus,  strength,  postyield  strain,  and  total  energy  were  calcu¬ 
lated. 

lateral,  antero-medial,  and  posterior  regions.  From  each  of 
these  regions,  one  to  three  prismatic  beams  were  cut  using 
a  diamond-coated  metallurgical  saw  (Isomet;  Buehler, 
Lake  Bluff,  IL,  USA).  The  beams  were  machined  to  regular 
test  samples  using  an  automated  CNC  milling  machine  un¬ 
der  constant  irrigation  (Modela  MDX-20;  Roland  DGA, 
Irvine,  CA,  USA).  Sample  width  (circumferential  direc¬ 
tion)  was  machined  to  5  mm  and  length  (longitudinal  di¬ 
rection)  was  machined  to  55  mm  for  all  samples.  Sample 
height  (radial  direction)  was  2.5  mm  except  for  four  tibias 
with  thin  cortices,  which  were  machined  to  2.2  mm.  A  total 
of  eight  samples  were  generated  from  each  tibia  and  ran¬ 
domly  distributed  to  monotonic  ( n  —  4)  and  damage  accu¬ 
mulation  ( n  —  4)  test  groups.  All  samples  were  stored  at 
-40°C  in  gauze  saturated  with  PBS  with  added  calcium*18-1 
and  placed  individually  in  airtight  bags. 

Monotonic  failure  properties 

Tissue-level  mechanical  properties  were  assessed  by 
loading  four  cortical  bone  samples  from  each  tibia  to  failure 
in  four-point  bending  at  0.05  mm/s  (Fig.  2A)  using  a  servo- 
hydraulic  materials  testing  system  (Instron  model  8872;  In- 
stron.  Canton,  MA,  USA).  Specimens  were  submerged  in  a 
PBS  solution  with  added  calcium*18-*  and  maintained  at 
37°C  throughout  all  tests.  Load  and  deflection  were  con¬ 
verted  to  stress  and  strain  using  the  following  equations, 
which  take  yielding  into  consideration*19*: 


a  =  2[2M  +  4>(dM/d4>)]/bh2  (2) 

s  =  h(j>/2a  =  !/2hA[(L  -  a)/(2a3/3  -  a2L  +  L3/3)]  (3) 

where  a  and  e  are  the  stress  and  strain  at  the  outer  surface 
of  the  beam,  M  =  applied  moment,  b  =  specimen  width,  h 
=  specimen  height,  a  =  V2  the  span  between  the  upper  two 
load  points  =  9  mm,  L  =  V2  the  span  between  the  two 
lower  load  points  =  21  mm,  =  angle  of  inclination  =  a/p, 
and  d/dcj>  is  the  derivative  with  respect  to  4>-  The  angle  of 
inclination  was  written  in  terms  of  the  measured  deflec¬ 
tion  (A)  by  estimating  the  curvature  (p)  using  standard 
beam  equations.  Mechanical  properties  were  calculated 
from  the  stress-strain  curves,  and  these  included  modulus, 
strength,  total  energy,  and  postyield  strain  (Fig.  2B).  Modu¬ 
lus  was  calculated  from  a  linear  regression  of  the  initial 
portion  of  the  stress-strain  curve.  Yield  was  determined 
using  the  0.2%  offset  method.  Postyield  strain  was  defined 
as  the  strain  at  failure  minus  the  strain  at  yield.  All  prop¬ 
erties  were  averaged  over  the  four  samples  tested  for  each 
tibia. 

Damage  accumulation  tests 

Tissue  damageability  was  assessed  using  a  protocol  de¬ 
signed  to  induce  and  accumulate  cracks  in  cortical  bone 
specimens.  The  accumulation  of  damage  leads  to  measur¬ 
able  degradation  of  mechanical  properties.*20-1  Therefore, 
the  degradation  of  mechanical  properties  can  be  used  as  an 
index  of  matrix  damage.  Four  cortical  bone  samples  from 
each  tibia  were  subjected  to  a  fifteen  cycle  damage  accu¬ 
mulation  protocol  (Fig.  3A)  similar  to  that  described  pre¬ 
viously.*21*  For  this  protocol,  “diagnostic”  cycles  (1,  3,  5,  7, 
9, 11, 13,  and  15)  were  interposed  between  “damage”  cycles 
(2,  4,  6,  8,  10,  12,  and  14).  For  the  diagnostic  cycles,  the 
specimens  were  loaded  in  four-point  bending  at  0.5  mm/s  to 
50%  of  the  average  displacement  at  yield  (determined  from 
the  monotonic  tests),  held  for  60  s,  and  unloaded  at  0.5 
mm/s.  Preliminary  studies  indicated  that  this  load  level  pro¬ 
vided  information  on  tissue-level  mechanical  properties 
without  inducing  additional  damage.  For  the  damage 
cycles,  the  specimens  were  loaded  at  0.5  mm/s  to  50%,  75%, 
100%,  125%,  150%,  175%,  and  200%  of  displacement  at 
yield,  respectively,  held  for  60  s,  and  unloaded  at  0.5  mm/s. 
A  5-minute  recovery  period  followed  each  damage  cycle. 
Displacement  at  yield  was  used  as  a  reference  in  the  dam¬ 
age  cycles  because  this  parameter  showed  little  variation 
among  the  test  samples  when  subjected  to  monotonic  four- 
point  bending.  The  displacement  at  yield  was  1.0  mm  for 
the  samples  with  a  height  of  2.5  mm  and  1.07  mm  for  the 
samples  with  a  height  of  2.2  mm. 

A  mechanical  measure  of  the  amount  of  damage  that 
accumulated  within  the  test  sample  was  quantified  from  the 
magnitude  of  stiffness  degradation.  For  each  diagnostic 
cycle,  stiffness  was  calculated  from  a  linear  regression  of  the 
initial  portion  of  the  load-deformation  curve.  Specimen 
stiffness  decreased  nonuniformly  with  each  cycle  revealing 
increasing  amounts  of  damage  induced  within  each  cycle 
and  an  overall  damage  accumulation  by  the  end  of  the  pro¬ 
tocol  (Fig.  3B).  At  the  end  of  the  test  sequence,  the  overall 
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A. 


Damage  Parameter:  D  =  1  -  S15  /  S0  D  =  0  (no  damage) 

D  =  1  (fracture) 


B. 


DAMAGE  CYCLE 

FIG.  3.  (A)  A  15-cycle  loading  protocol  was  used  to  induce  dam¬ 

age  within  machined  cortical  bone  specimens  and  to  measure  re¬ 
sultant  stiffness  degradation.  Damage  was  induced  during  cycles 
2, 4,  6,  8, 10, 12,  and  14  by  conducting  relaxation  tests  at  increasing 
levels  of  applied  displacement  (expressed  as  a  percentage  of  the 
displacement  at  yield).  Diagnostic  cycles  were  interposed  between 
damaging  cycles.  Each  damage  cycle  was  preceded  by  a  diagnostic 
cycle  at  50%  of  displacement  at  yield.  A  5-minute  recovery  period 
was  introduced  after  the  damage  cycles  to  relieve  residual  internal 
stresses.  (B)  The  change  in  stiffness  calculated  between  sequential 
diagnostic  tests  was  plotted  vs.  cycle  number  for  the  damage  tests. 
The  dashed  curve  represents  a  specimen  showing  little  stiffness 
degradation  (i.e.,  little  damage  accumulation).  The  solid  line  rep¬ 
resents  a  specimen  showing  large  stiffness  degradation  (i.e.,  more 
damage  accumulation). 


damage  parameter,  D ,  was  calculated  by  comparing  the 
stiffness  of  the  first  and  last  diagnostic  tests  such  that: 

D  =  1  -  Sls/S0,  (4) 

where  S15  is  the  stiffness  of  the  last  diagnostic  cycle  and  S0 
is  the  average  stiffness  of  the  first  two  diagnostic  cycles  (Sj, 
S3)  and  the  first  damage  cycle  ( S2 )■ 

Statistical  analysis 

All  data  were  regressed  against  age  using  linear  regres¬ 
sion  analysis  to  identify  the  properties  that  varied  signifi¬ 
cantly  with  age  (GraphPad  Prism,  San  Diego,  CA,  USA). 
To  determine  whether  bone  morphology  was  related  to  tis¬ 
sue  level  material  properties,  partial  correlation  coefficients 
were  determined  between  each  geometrical  parameter 
(e.g.,  IA P,  /ML,  J,  S)  and  each  tissue  level  mechanical  prop¬ 
erty  (modulus,  strength,  total  energy,  postyield  strain,  dam- 
ageability)  while  taking  age  into  consideration  (Minitab, 
State  College,  PA,  USA).(22) 

RESULTS 

The  sample  population  showed  broad  ranges  of  body 
size,  body  stature,  and  bone  morphology  values  (Table  1). 
Modulus  and  strength  showed  little  variation  among  indi¬ 
viduals  (CV  =  9.73%  and  4.62%,  respectively).  However, 


Table  1.  Variation  in  Properties  Among  Young  Adult 
Male  Tibias 


Property  Mean  ±  SD  Range 


Age  (years) 

Body  weight  (kg) 

Body  height  (cm) 

Body  mass  index  (kg/m2) 

Tibia  length  (cm) 

Cortical  area  (mm2) 

AP  width  (mm) 

ML  width  (mm) 

AP  moment  of  inertia  (mm4) 
ML  moment  of  inertia  (mm4) 
Polar  moment  of  inertia  (mm4) 
AP  section  modulus  (mm3) 

ML  section  modulus  (mm3) 

AP  slenderness  (l/mm2/kg) 

ML  slenderness  (l/mm2/kg) 
Modulus  (GPa) 

Strength  (MPa) 

Total  energy  (MPa) 

Post-yield  strain 
Damage  parameter  (D) 


32.9  ±  10.4 

17-46 

83.8  ±  25.1 

57.2-158.8 

177.7  ±  4.3 

170.2-182.9 

26.7  ±  8.0 

17.1-51.7 

38.1  ±  1.9 

34.4-40.4 

355.9  ±  55.2 

268.3-511.1 

31.2  ±2.5 

28.4-36.9 

24.3  ±  2.3 

19.9-30.6 

34,390  ±  10,149 

19,466-60,698 

17,250  ±  5,945 

8,246-34,734 

51,640  ±15,886 

27,713-95,432 

3,279  ±819 

1,925-5,172 

4,188  ±  907 

2,785-6,237 

9.9  ±  2.0 

6.5-13.4 

7.7  ±1.5 

5.4-10.2 

17.0  ±1.7 

13.4-19.0 

130.7  ±  6.1 

120.8-144.5 

3.3  ±  0.9 

2.3-5. 6 

0.026  ±  0.006 

0.016-0.039 

0.165  ±  0.038 

0.103-0.253 

AP,  antero-posterior  direction;  ML,  medial-lateral  direction. 


postyield  strain  (CV  =  24.0%),  total  energy  (CV  = 
26.4%),  and  the  damage  parameter  (CV  =  23.0%)  all 
showed  large  variability  among  the  samples.  Morphological 
measures  such  as  AP  width,  section  modulus,  and  the  polar 
moment  of  inertia,  J  (Fig.  4),  increased  linearly  with  body 
weight  (R2  —  0.59,  p  <  0.003)  and  body  mass  index  (BMI; 
R2  —  0.57,  p  <  0.004),  but  were  independent  of  body  height 
(R2  =  0.01,  p  <  0.7).  These  relationships  did  not  change 
when  the  body  weight  values  were  corrected  for  age  (data 
not  shown).  Body  height  was  uncorrelated  with  body 
weight  (R2  =  0.01,  p  <  0.8),  indicating  that  the  sample 
population  consisted  of  individuals  with  similar  heights  but 
widely  varying  body  weights. 

Significant  age-related  changes  were  observed  for  the  tis¬ 
sue-level  mechanical  properties  and  the  size  of  the  tibia.  A 
significant,  positive  correlation  was  observed  between  tibia 
slenderness  in  the  AP  (R2  =  0.31,  p  <  0.02)  and  ML  (R2  = 
0.24,  p  <  0.05)  directions  and  age.  However,  /AP,  /ML,  and 
J  did  not  vary  with  age,  suggesting  that  the  variation  in 
slenderness  with  age  was  due  largely  to  higher  body  weight 
and  BMI  (R2  =  0.29-0.32,  p  <  0.03)  values  for  the  older 
individuals.  Although  tissue  modulus  did  not  vary  signifi¬ 
cantly  with  age,  tissue  strength  (R2  =  0.53,  p  <  0.001),  post¬ 
yield  strain  (R2  —  0.44,  p  <  0.004),  and  total  energy  (R2  = 
0.32,  p  <  0.002)  were  significantly  lower  for  the  older  indi¬ 
viduals.  Furthermore,  a  significant,  negative  correlation  was 
observed  between  the  damage  parameter  and  age  (R2  — 
0.41,  p  <  0.006).  This  data  suggested  that,  whereas  the  tibia 
became  more  slender  relative  to  body  size  with  age,  the 
cortical  tissue  became  progressively  less  strong  and  less 
ductile  (i.e.,  more  brittle)  with  age. 

The  correlation  analysis  showed  significant  correlations 
between  tibial  morphology  and  the  mechanical  properties 
that  characterized  tissue  brittleness  and  damageability 
(Table  2).  The  relationships  among  tissue-level  mechanical 
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FIG.  4.  The  average  cross-sectional  polar  moment  of  inertia  of 
the  tibia  from  males  increased  with  (A)  body  weight  and  (B)  body 
mass  index  but  not  with  (C)  body  height.  Removing  the  outlier 
changed  the  R2  values  to  0.24  (p  <  0.05)  for  A  and  0.20  (p  <  0.09) 
for  B. 


properties  and  cross-sectional  morphology  were  linear.  Post¬ 
yield  strain  and  total  energy  increased  significantly  with 
AP  width  (Figs.  5 A  and  5B).  Modulus  decreased  with  7AP 
(p  <  0.07),  J  (p  <  0.08),  AP  section  modulus  (p  <  0.05),  and 
ML  section  modulus  (Figs.  5C-5F).  Tissue  damageability 
increased  with  tibia  slenderness  in  the  AP  (p  <  0.05;  Fig.  6) 
and  ML  (p  <  0.09)  directions.  These  correlations,  which 
were  independent  of  age,  indicated  that  a  narrower  bone 


was  comprised  of  tissue  that  failed  in  a  more  brittle  man¬ 
ner  and  accumulated  more  damage. 

DISCUSSION 

The  results  of  this  study  revealed  that  the  tissue-level 
mechanical  properties  of  cortical  bone  varied  with  the  size 
of  the  tibia.  Positive  correlations  were  observed  between 
measures  of  bone  size  (AP  width)  and  measures  of  tissue 
ductility  (postyield  strain,  total  energy),  and  negative  cor¬ 
relations  were  observed  between  bone  size  (moment  of  in¬ 
ertia,  section  modulus)  and  tissue  modulus.  Many  of  these 
correlations  were  significant.  The  lack  of  significant  corre¬ 
lation  with  all  measures  of  bone  size  can  be  attributed 
largely  to  the  complex  shape  of  the  tibia.  The  tibia  has  a 
triangular  cross-section  and,  consequently,  measures  of 
width  correlated  significantly  with  mechanically  relevant 
traits  like  cortical  area  and  moment  of  inertia  but  explained 
only  50-80%  of  the  variability  in  these  measures  (data  not 
shown).  These  correlations  would  be  greater  if  the  cross- 
section  had  a  circular  shape.  The  variability  in  these  corre¬ 
lations  was  sufficiently  large  that  neither  the  linear  (width) 
traits  nor  the  integrated  traits  like  area  and  moment  of 
inertia  correlated  significantly  with  a  particular  tissue -level 
mechanical  property  simultaneously.  Nevertheless,  the  data 
indicated  that  bones  with  smaller  width  were  comprised  of 
stiffer  and  less  ductile  (i.e.,  more  brittle)  material  compared 
with  larger,  more  robust  bones.  The  correlation  between 
tissue  ductility  and  bone  size  may  help  explain  why  male 
military  recruits^1,2,11-1  and  male  athletes(12)  with  narrow 
bones  show  a  higher  incidence  of  stress  fractures  compared 
with  individuals  with  wide  bones. 

The  development  of  the  slenderness  index(17)  was  for  a 
“normal”  range  in  height  and  weight  and  is  probably  not 
useful  beyond  this  range.  Flowever,  the  morphological 
variation  observed  in  our  sample  population  was  consistent 
with  that  reported  for  military  recruits(1,2)  and  runners/12-1 
and  height  and  weight  were  consistent  with  recent  national 
averages. (23)  As  expected,  bone  size  varied  with  body 
weight/17-1  but  did  not  vary  with  height  (Fig.  4)/24>  Thus, 
narrow  bones  came  from  less  heavy  individuals  who  were  of 
similar  height  as  those  with  wide  tibias.  Weight  varied  more 
than  height  for  our  sample  population  similar  to  that  ob¬ 
served  for  the  aged-matched  national  data.  Furthermore, 
the  variability  in  weight,  specifically  inclusion  of  one  outlier 
(Fig.  4),  did  not  affect  the  results  (i.e.,  the  heaviest  person 
did  not  have  an  unusual  slenderness  value).  Thus,  the  bones 
used  in  this  study  seem  to  be  an  appropriate  size  relative  to 
body  type. 

The  variation  in  long  bone  slenderness  has  been  attrib¬ 
uted  to  genetic  and  environmental  factors  influencing 
growth  and  development^  and  has  been  implicated  as  a 
risk  factor  for  osteoporotic  fracture/26)  To  be  relevant  for 
military  recruits,  the  sample  population  should  have  ranged 
in  age  between  18  and  25  years.  However,  for  the  age  range 
in  this  study,  the  tissue-level  mechanical  properties  varied 
linearly  with  age  and  were  easily  corrected  using  a  linear 
regression  method/22-1  Consequently,  the  correlation  analy¬ 
sis  presented  here  provides  relevant  insight  into  the  rela¬ 
tionship  observed  between  bone  size  and  stress  fracture  risk 
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Table  2.  Partial  Correlation  Coefficients  Taking  Age  Into  Consideration 


Modulus 

Strength 

Total  energy 

PY  strain 

Damageability 

Cortical  area 

-0.24 

0.03 

0.37 

0.40 

-0.23 

(0.35) 

(0.90) 

(0.14) 

(0.11) 

(0.37) 

AP  width 

-0.09 

-0.03 

0.57 

0.70 

-0.16 

(0.74) 

(0.90) 

(0.02) 

(0.002) 

(0.55) 

ML  width 

-0.32 

-0.22 

0.34 

0.41 

-0.19 

(0.21) 

(0.39) 

(0.18) 

(0.11) 

(0.45) 

AP  moment  of  inertia 

-0.45 

-0.10 

0.25 

0.32 

-0.27 

(0.07) 

(0.70) 

(0.34) 

(0.21) 

(0.29) 

ML  moment  of  inertia 

-0.39 

-0.18 

0.22 

0.28 

-0.25 

(0.12) 

(0.50) 

(0.39) 

(0.28) 

(0.32) 

Polar  moment  of  inertia 

-0.43 

-0.13 

0.24 

0.31 

-0.27 

(0.08) 

(0.61) 

(0.35) 

(0.23) 

(0.30) 

AP  section  modulus 

-0.50 

-0.13 

0.10 

0.14 

-0.30 

(0.04) 

(0.62) 

(0.70) 

(0.60) 

(0.25) 

ML  section  modulus 

-0.47 

-0.05 

0.18 

0.23 

-0.33 

(0.06) 

(0.84) 

(0.49) 

(0.38) 

(0.20) 

Ap  slenderness 

0.36 

0.10 

0.09 

0.03 

0.48 

(0.15) 

(0.69) 

(0.72) 

(0.92) 

(0.05) 

ML  slenderness 

0.23 

-0.01 

0.04 

-0.05 

0.41 

(0.39) 

(0.96) 

(0.89) 

(0.84) 

(0.10) 

Pearson  correlation  coefficients 

are  shown  with  p  values 

in  parentheses.  Significant  correlations  are 

shown  in  bold.  Abbreviations 

are  as  shown  in 

Table  1. 


for  young  adult  males.  Further  studies  are  needed  to  deter¬ 
mine  if  this  relationship  holds  over  a  wider  (older)  age 
range. 

The  data  provide  a  new  paradigm  that  may  explain  how 
variation  in  bone  slenderness  contributes  to  stress  fracture 
risk.  Individuals  with  narrow  tibias  were  previously  thought 
to  show  increased  fatigue  damage  during  intense  training 
because  the  smaller  bone  size  would  lead  to  an  overload 
situation  (i.e.,  higher  tissue  level  stresses).' (1,2,12:i  This  inter¬ 
pretation  was  based  on  the  assumption  that  tissue  mechani¬ 
cal  properties  did  not  vary  among  individuals.  However,  the 
current  results  indicated  that  tissue-level  mechanical  prop¬ 
erties  do  vary  among  individuals.  Specifically,  the  data  sug¬ 
gest  that  there  are  at  least  two  important  tissue-level  me¬ 
chanical  property  variations  that  need  to  be  considered  to 
understand  why  bone  size  is  a  risk  factor  for  stress  fractures. 
Narrower  tibias  were  comprised  of  tissue  that  was  more 
brittle  (low  total  energy)  and  was  prone  to  accumulate 
more  damage  compared  with  tissue  from  wider  tibia.  Hav¬ 
ing  tissue  that  is  more  or  less  damageable  may  be  inconse¬ 
quential  during  day-to-day  activities.  However,  tissue-level 
mechanical  properties  like  total  energy  and  ductility  be¬ 
come  particularly  important  in  defining  the  response  of 
bone  to  an  extreme  loading  condition,  such  as  that  expected 
during  military  training  or  during  a  fall.  Total  energy  de¬ 
fines  the  amount  of  energy  required  to  break  a  bone  (im¬ 
portant  during  a  fall)  and  ductility  and  damageability  define 
the  amount  of  damage  accumulated  under  overload  or  re¬ 
petitive  loading  (important  during  military  training).  Fur¬ 
thermore,  tissue  stresses  would  be  expected  to  remain 
higher  for  narrow  tibias  loaded  in  bending  or  torsion.  Mo¬ 
ment  of  inertia  is  related  to  the  external  diameter  raised  to 
the  fourth  power.  Because  whole  bone  stiffness  and 
strength  are  correlated  with  moment  of  inertia/17,27*  a  bone 
with  a  large  external  diameter  should  also  show  large  over¬ 


all  stiffness  and  strength  values.  However,  the  -30%  varia¬ 
tion  in  tissue  modulus  (Table  1)  did  not  fully  compensate 
for  the  -100%  variation  in  the  moment  of  inertia  or  the 
section  modulus  (Table  l)/27*  Thus,  in  situ  damage  accu¬ 
mulation  may  elicit  a  biological  response  (remodeling)  that, 
coupled  with  the  higher  tissue  stresses,  exacerbates  the  fa¬ 
tigue  process/13,28*  Consequently,  individuals  with  narrow 
tibia  may  be  at  higher  risk  of  stress  fractures  because  of 
higher  in  vivo  tissue  stresses  (overloading)  coupled  with 
tissue  that  is  more  prone  to  accumulating  damage. 

The  data  may  also  help  explain  why  age  is  another  risk 
factor  for  stress  fractures/7,29-1  Bone  strength,  postyield 
strain,  and  total  energy  decreased  over  the  17-  to  46-year-old 
age  range.  This  was  consistent  with  previous  studies130-32* 
and  indicated  that  cortical  bone  becomes  less  ductile  (i.e., 
more  brittle)  and  weaker  with  age  and  that  these  changes 
began  early  in  life.  This  age-related  decline  in  strength  and 
ductility  is  thought  to  be  a  result  of  increased  mineraliza¬ 
tion  and  remodeling/33,34-1  Thus,  even  in  the  young  adult 
age  range,  the  amount  of  damage  accumulated  under  vig¬ 
orous  loading  regimens  would  be  expected  to  increase  with 
age.  This  variation  in  tissue  ductility  may  increase  the  sus¬ 
ceptibility  of  stress  fracture  risk  for  recruits  that  enter  into 
military  training  at  an  older  age. 

These  results,  and  those  of  others,  indicated  that  not 
all  cortical  tissue  was  constructed  in  the  same  manner. 
The  mechanical  properties  of  cortical  tissue  vary  with 
age,130-32,35*  across  species/35*  among  bones  of  the  same 
individual/36*  and  among  anatomical  sites  within  the  same 
bone/37,38*  Here  we  showed  that  the  average  mechanical 
properties  of  cortical  tissue  also  varied  as  a  function  of  the 
overall  size  of  the  bone.  This  coupling  between  bone  mor¬ 
phology  and  tissue -level  mechanical  properties  has  been 
attributed  to  an  adaptive  response  of  bone/14,35,39*  In  this 
study,  smaller  tibia  bone  size  was  coupled  with  an  increase 
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FIG.  5.  (A)  Postyield  strain  and  (B)  total 
energy  correlated  with  AP  width.  Modulus 
decreased  with  (C)  /AP,  (D)  J,  (E)  AP  sec¬ 
tion  modulus,  and  (F)  ML  section  modulus. 
Data  were  age-corrected  based  on  a  linear 
regression  method/22* 


Robust  Gracile 
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FIG.  6.  Damageability  correlated  with  AP  slenderness  suggest¬ 
ing  that  tibias  that  were  more  slender  relative  to  body  size  and 
stature  were  comprised  of  tissue  that  accumulated  more  damage. 
Data  were  age-corrected  based  on  a  linear  regression  method.'22-1 

in  tissue  modulus.  The  goal  of  this  adaptive  response  is  to 
ensure  that  morphology  and  quality  together  meet  me¬ 
chanical  demands.  This  coupling  was  observed  when  com¬ 
paring  bones  subjected  to  widely  varying  mechanical  de¬ 
mands  from  different  species^35*  and  has  also  been  used  to 
explain  the  maturation  of  bone  during  growth/40-43*  Our 


current  results  suggested  that  this  coupling  might  also  exist 
for  a  particular  bone  (tibia)  within  the  same  species  (hu¬ 
man).  Additional  studies  are  needed  to  determine  if  similar 
relationships  between  morphology  and  quality  exist  for 
other  long  bones  (femur,  humerus,  radius). 

The  relationship  between  morphology  and  tissue -level 
mechanical  properties  observed  in  the  human  skeleton  was 
consistent  with  that  observed  for  the  mouse  skeleton/14*  In 
both  the  mouse  and  human  skeletons,  genetic  heterogene¬ 
ity  leads  to  variability  in  adult  bone  morphology  and  tissue 
level  mechanical  properties.  A  comparison  of  femurs  from 
A/J  and  C57/BL6  (B6)  inbred  strains  showed  that  A/J  fe¬ 
murs  were  more  slender  than  B6  as  a  result  of  the  two 
strains  having  similar  bone  lengths,  but  A/J  having  a  sig¬ 
nificantly  smaller  cross-sectional  size  and  shape/14*  Despite 
the  difference  in  bone  size,  the  two  strains  showed  similar 
whole  bone  stiffness  values.  The  variability  in  bone  slen¬ 
derness  was  inversely  related  to  mineral  content,  suggesting 
that  bone  morphology  and  mineral  content  were  coordi- 
nately  regulated  so  whole  bone  stiffness  appropriately 
matched  the  mechanical  demands  imposed  by  weight  bear¬ 
ing.  However,  as  a  result  of  regulating  mineral  content  to 
match  bone  size,  A/J  femurs  failed  in  a  brittle  manner  and 
showed  poor  fatigue  properties.  In  the  human  skeleton, 
smaller  bones  were  stiffer  and  less  ductile.  Thus,  a  recipro¬ 
cal  relationship  was  observed  between  bone  stiffness  and 
ductility  for  both  skeleton  systems.  This  reciprocal  relation¬ 
ship  has  been  extensively  reported  for  cortical  bone/15*  and 
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it  is  thought  to  be  a  result  of  the  nature  of  the  compositional 
and  structural  factors  that  can  be  modulated  on  a  biological 
level.'44^47 'Although  variation  in  mineral  content  may  have 
explained  the  differences  in  brittleness  for  the  mouse  skel¬ 
eton,  we  expect  that  the  human  skeleton  will  be  more  com¬ 
plex  and  that  the  variation  in  tissue -level  mechanical  prop¬ 
erties  will  be  a  consequence  of  variable  composition 
(mineral,  collagen,  water)  as  well  as  microarchitecture  (la¬ 
mellae,  osteon  size,  porosity). 

The  calculated  bending  modulus  and  strength  values, 
which  were  determined  from  machined  bone  samples  and 
were  thus  quantified  in  a  manner  that  was  independent  of 
bone  size,  were  consistent  with  bone  tensile  properties/30' 
as  expected.  Test  samples  were  randomly  selected  to  obtain 
representative  mean  values  for  each  tibia  and  the  variation 
in  mechanical  properties  within  each  tibia  was  similar  to  the 
variability  observed  across  tibias.  Thus,  we  believe  that  the 
mean  values  reported  here  represent  the  generalized  tissue- 
level  mechanical  behavior  for  each  tibia. 

Compared  with  back-calculating  tissue-level  mechanical 
properties  from  whole  bone  failure  tests,  the  current 
method  of  measuring  tissue-level  mechanical  properties  di¬ 
rectly  from  machined  samples  provided  a  broader  range  of 
mechanical  properties  that  were  needed  to  better  under¬ 
stand  why  bone  size  is  a  risk  factor  for  stress  fractures.  The 
mechanical  properties  included  measures  of  ductility  (i.e., 
postyield  strain,  total  energy)  as  well  as  an  independent 
measure  of  damageability  (i.e.,  the  damage  parameter). 
These  properties  were  chosen  because  they  were  relevant 
for  understanding  the  material  response  of  bones  subjected 
to  the  vigorous,  repetitive  loading  associated  with  military 
training  and  running.  Postyield  strain  and  total  energy  rep¬ 
resent  measures  of  tissue  ductility  and  were  assessed  to 
discriminate  between  ductile  and  brittle  failure  modes.  Ma¬ 
terials  that  fail  in  a  brittle  manner  show  low  postyield  strain 
and  total  energy  values.  Variation  in  the  ductility  of  cortical 
bone  arises  from  differences  in  the  initiation,  accumulation, 
propagation,  and  coalescence  of  damage  in  the  form  of  mi¬ 
crocracks/48'49'  Variation  in  the  damage  parameter  re¬ 
flected  differences  in  the  amount  of  damage  accumulated 
within  the  tissue  and/or  differences  in  the  way  damage  de¬ 
graded  tissue  stiffness.  The  damage  parameter  correlated 
negatively  with  postyield  strain  and  total  energy  (R2  — 
0.22-0.25 ,p  <  0.05)  indicating  that  cortical  tissue  that  failed 
in  a  brittle  manner  also  tended  to  have  higher  tissue  dam¬ 
ageability  or,  accumulate  more  damage.  Although  the  ex 
vivo  bending  tests  do  not  necessarily  reflect  the  in  vivo 
loads  imposed  on  the  tibia,1 (17’5°-52)  the  bending  loads  were 
expected  to  induce  a  combination  of  tensile,  compressive, 
and  shear  damage'53'  that  may  be  sufficiently  complex  to 
represent  a  generalized  variation  in  bone  quality  among 
human  tibias. 

The  results  of  this  study  provide  new  insight  into  why 
bone  size  is  a  risk  factor  for  stress  fractures.  Stress  fractures 
are  believed  to  be  a  consequence  of  excess  damage  accu¬ 
mulation  following  intense,  repetitive  activities.  Biological 
processes  that  attempt  to  repair  the  damage  may  further 
weaken  the  tissue  because  the  increased  resorption  results 
in  increased  tissue  porosity.'54'  However,  the  actual  contri¬ 
bution  of  biological  repair  processes  to  stress  fracture  risk 


remains  unclear.'55'  Damage,  in  the  form  of  microcracks,  is 
the  expected  sequelae  of  repetitive  loading  following  nor¬ 
mal,  daily  activities.'56'  Intense  loading  conditions,  such  as 
those  associated  with  military  training  and  long  distance 
running,  are  expected  to  further  increase  in  situ  damage 
accumulation  and  degrade  tissue-level  mechanical  proper¬ 
ties.'13'  Therefore,  under  extreme  loading  conditions  (e.g., 
military  training),  variation  in  bone  quality,  specifically  tis¬ 
sue  damageability,  may  be  a  contributing  factor  to  the  in¬ 
creased  risk  of  stress  fracture  for  individuals  with  more 
slender  bones.  The  current  data  suggested  that  bone  mor¬ 
phology  could  be  used  as  a  predictor  of  tissue  fragility  and 
stress  fracture  risk  in  the  absence  of  available  noninvasive 
imaging  techniques  that  accurately  measure  bone  damage- 
ability. 
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Genetic  Variation  in  Bone  Growth  Patterns  Defines  Adult  Mouse 

Bone  Fragility 
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ABSTRACT:  Femoral  morphology  and  composition  were  determined  for  three  inbred  mouse  strains  between 
ages  E18.5  and  1  year.  Genotype-specific  variation  in  postnatal,  pubertal,  and  postpubertal  growth  patterns 
and  mineral  accrual  explained  differences  in  adult  bone  trait  combinations  and  thus  bone  fragility. 

Introduction:  Fracture  risk  is  strongly  regulated  by  genetic  factors.  However,  this  regulation  is  generally 
considered  complex  and  polygenic.  Therefore,  the  development  of  effective  genetic-based  diagnostic  and 
treatment  tools  hinges  on  understanding  how  multiple  genes  and  multiple  cell  types  interact  to  create  me¬ 
chanically  functional  structures.  The  goal  of  this  study  was  to  connect  variability  in  whole  bone  mechanical 
function,  including  measures  of  fragility,  to  variability  in  the  biological  processes  underlying  skeletal  devel¬ 
opment.  We  accomplished  this  by  testing  for  variation  in  bone  morphology  and  composition  among  three 
inbred  mouse  strains  from  E18.5  to  1  year  of  age. 

Materials  and  Methods:  Mid-diaphyseal  cross-sectional  areas,  diameters,  moments  of  inertia,  and  ash  content 
were  determined  for  three  strains  of  mice  with  widely  differing  adult  whole  bone  femoral  mechanical  prop¬ 
erties  (A/J,  C57BL/6J,  and  C3H/HeJ)  at  E18.5  and  postnatal  days  1,  7, 14,  28,  56, 112, 182,  and  365  (n  =  5—15 
mice/strain/age). 

Results:  Significant  differences  in  the  magnitude  and  rate  of  change  in  morphological  and  compositional  bone 
traits  were  observed  among  the  three  strains  at  each  phase  of  growth,  including  prenatal,  postnatal,  pubertal, 
and  adult  ages.  These  genotype-specific  variations  in  growth  patterns  mathematically  determined  how  varia¬ 
tion  in  adult  bone  trait  combinations  and  mechanical  properties  arose.  Furthermore,  six  bone  traits  were 
identified  that  characterize  phenotypic  variability  in  femoral  growth.  These  include  (1)  bone  size  and  shape  at 
postnatal  day  1,  (2)  periosteal  and  (3)  endosteal  expansion  during  early  growth,  (4)  periosteal  expansion  and 
(5)  endosteal  contraction  in  later  growth,  and  (6)  ash  content.  These  results  show  that  genetic  variability  in 
adult  bone  traits  arises  from  variation  in  biological  processes  at  each  phase  of  growth. 

Conclusions:  Inbred  mice  achieve  different  combinations  of  adult  bone  traits  through  genotype-specific  regu¬ 
lation  of  bone  surface  activity,  growth  patterns,  and  whole  bone  mineral  accrual  throughout  femoral  devel¬ 
opment.  This  study  provides  a  systematic  approach,  which  can  be  applied  to  the  human  skeleton,  to  uncover 
genetic  control  mechanisms  influencing  bone  fragility. 
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INTRODUCTION 

Genetic  factors  are  generally  acknowledged  to  play 
an  important  role  in  the  determination  of  skeletal  fra¬ 
gility.  Clinically,  genetic  variability  in  BMD,(1)  bone  qual¬ 
ity, (2)  adult  bone  morphology,®  and  the  kinetics  of  bone 
loss<4  5>  all  contribute  to  the  variation  in  fracture  risk/61 
Therefore,  genetic-based  screening  tools  that  identify  an 
individual’s  skeletal  growth  potential  and/or  potential  for 
bone  loss  could  be  used  for  the  early  diagnosis  of  osteopo¬ 
rotic  fracture  risk  and  for  tailored  treatment  regimens 
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thereafter/71  However,  these  genetic  factors  are  not  com¬ 
pletely  understood.  Additionally,  the  processes  of  bone 
growth  and  bone  loss  are  genetically  complex,  such  that  the 
modeling  and  remodeling  of  complex  structures  involves 
the  coordinated  expression  of  many  genes  by  multiple  cell 
types  (osteoblasts,  osteoclasts,  and  osteocytes)/7’  There¬ 
fore,  a  primary  challenge  for  developing  effective  genetic- 
based  diagnostic  tools  for  fracture  risk  will  be  to  understand 
how  multiple  genes,  in  multiple  cell  types,  interact  through¬ 
out  life  to  create  and  maintain  structures  that  are  mechani¬ 
cally  and  biologically  functional.® 

To  study  the  functional  connections  between  genotype 
and  skeletal  fragility,  we  adopted  a  top-down  approach  that 
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reveals  the  functional  connections  between  genotype  and 
skeletal  fragility/9-’  Our  approach  uses  long  bones  from  in- 
bred  mouse  strains  to  relate  genetic  variability  in  whole 
bone  mechanical  properties  to  the  underlying  biological 
processes  regulating  bone  growth  and  loss  in  a  systematic 
and  hierarchical  manner.  Previous  research  using  inbred 
mice  has  shown  genotype-specific  variability  in  skeletal 
morphology,1 (10-12’  BMD,(13)  mechanical  properties/14,15’ 
bone  formation  rates/16,17’  and  responsiveness  to  mechani¬ 
cal  stimulation/18,19’  Previously,  we  showed  that  individual 
inbred  mouse  strains  exhibited  a  particular  combination  of 
cortical  area  (Ct.Ar),  polar  moment  of  inertia  (JD),  and 
mineral  content  (ash  content)  that  together  explained  66- 
88%  of  the  genetic  variability  in  adult  whole  bone  mechani¬ 
cal  properties/9’  These  results  indicated  that  the  combina¬ 
tion  of  transverse  bone  size,  shape,  and  material  quality  was 
responsible  for  the  particular  repertoire  of  stiffness, 
strength,  and  toughness  that  characterized  each  strain.  The 
next  step  in  our  systematic  approach  was  to  determine  how 
genetic  variability  in  these  three  physical  bone  traits  (Ct.Ar, 
J0,  and  ash  content)  arises  during  the  growth  process. 

Bone  traits  are  influenced  by  both  genetic  and  environ¬ 
mental  factors  during  growth,  and  specific  combinations  of 
bone  traits  are  established  to  satisfy  the  mechanical  de¬ 
mands  associated  with  weight  bearing/20’  Long  bone  di- 
aphyses  are  hollow  cylinders  whose  size  and  shape  can  be 
mathematically  determined  by  the  relative  amounts  of  bone 
apposition  and  resorption  on  the  periosteal  and  endosteal 
surfaces/21-24’  We  hypothesize  that  genetic  variability  in 
the  rates  of  periosteal  and  endosteal  expansion  throughout 
life  (i.e.,  growth  patterns)  defines  genotype-specific  combi¬ 
nations  of  adult  bone  traits.  To  test  this  hypothesis,  we 
compared  changes  in  femoral  morphology  (Ct.Ar  and  JD) 
and  composition  (ash  content)  across  development  (from 
embryonic  day  18.5  [E18.5]  to  1  year).  Transverse  growth 
was  studied  because  these  traits  contribute  predominantly 
to  the  resistance  of  in  vivo  loads/25’  Three  genetically  dis¬ 
tinct  inbred  mouse  strains,  A/J  (A),  C57BL/6J  (B6),  and 
C3H/HeJ  (C3H),  with  widely  differing  adult  bone  trait 
magnitudes,  were  examined  to  better  understand  the  ex¬ 
pansive  growth  of  bone  surfaces  and  how  this  growth  leads 
to  unique  combinations  of  adult  morphological  (Ct.Ar, 
Tt.Ar,  Ma.Ar,  and  JD)  and  compositional  (ash  content)  trait 
values/9’  We  show  that  variability  in  skeletal  growth  pat¬ 
terns  reveals  the  biological  and  genetic  control  mechanisms 
that  define  adult  physical  bone  traits  and,  consequently, 
bone  fragility. 

MATERIALS  AND  METHODS 

Animals 

Female  (n  =  20/strain)  and  male  ( n  =  10/strain)  A/J  (A), 
C57BL/6J  (B6),  and  C3H/HeJ  (C3H)  mice  were  purchased 
from  Jackson  Laboratory  (Bar  Harbor,  ME,  USA)  at  6-8 
weeks  of  age  and  used  to  establish  breeding  colonies.  Mice 
were  fed  a  standard  mouse  chow  (Purina  Laboratory  Chow 
5001;  Purina  Mills)  and  water  ad  libitum  and  kept  on  a  12-h 
light:dark  cycle.  Matings  were  monitored,  and  after  4  weeks 
of  weaning,  the  offspring  were  separated  by  sex  and  housed 


with  four  to  five  mice  per  cage.  Female  offspring  were 
killed  at  selected  time -points  for  analysis:  these  included 
E18.5,  postnatal  day  1, 7, 14, 28,  56, 112, 182,  and  365  (Table 
1).  For  the  E18.5  age  group,  pregnant  females  were  identi¬ 
fied  by  the  first  appearance  of  a  vaginal  plug  (designated  as 
day  E0.5).  The  pups  were  extracted  18  days  later  through 
caesarian  section  and  killed  by  decapitation.  All  mice  1  day 
of  age  and  older  were  killed  by  carbon  dioxide  asphyxiation 
and  decapitation.  The  Mount  Sinai  Institutional  Animal 
Care  and  Use  Committee  approved  all  procedures  for  the 
treatment  of  mice. 

Morphological  traits 

Right  femurs  were  harvested  and  embedded  in  poly¬ 
methylmethacrylate  as  described  previously/14’  Transverse 
sections  of  the  femoral  diaphyses  (150  p,m  in  thickness) 
were  obtained  using  a  low-speed  diamond-coated  wafering 
saw  (Buehler,  Lake  Bluff,  IL,  USA)  from  a  site  immedi¬ 
ately  distal  to  the  third  trochanter.  The  sections  were  af¬ 
fixed  to  acrylic  slides,  polished  to  1-pm  finish,  and  stained 
using  von  Kossa’s  method/26’  The  sections  were  imaged 
with  a  digital  camera  (Sony  Exwave  HAD,  3CCD  Camera; 
Sony)  attached  to  a  visible  light  microscope  (Zeiss  Ax- 
ioplan2;  Zeiss).  The  pixel  resolution  of  this  optical  imaging 
system  was  2.1  p,m.  An  advanced  image  analysis  software 
package  (IMAQ  Vision  Builder  6.0;  National  Instruments 
Corp.,  Austin,  TX,  USA)  was  used  to  quantify  the  cross- 
sectional  areas  (total  area  [Tt.Ar],  cortical  area  [Ct.Ar],  and 
marrow  area  [Ma.Ar]),  cortical  thickness  (Ct.Th),  and  the 
polar  moment  of  inertia  (J0).  Young  bones  (ages  E18.5 
through  2  weeks  of  age)  had  a  predominantly  woven  ap¬ 
pearance  with  considerable  cortical  porosity.  Therefore,  at 
these  time -points,  morphometric  measurements  included 
all  porosities  within  the  cortical  shell  and  were  considered 
apparent  traits.  Three  cross-sections  were  analyzed  per 
bone,  and  the  values  were  averaged.  To  establish  rates  of 
surface  growth,  the  time  normalized  change  in  Tt.Ar  and 
Ma.Ar  during  the  early  phase  of  growth  (1-28  days  of  age) 
and  the  later  phase  of  growth  (28-112  days  of  age)  were 
obtained  from  the  following  equations  for  all  strains: 

ATt.Ar]_28  =  [Tt.Ar28  -  avg(Tt.Ar1)]/28  days 

(1) 

ATt.Ar28_112  =  [Tt.Ar112  -  avg(Tt.Ar28)]/84  days 

(2) 

AMa.Ar]_28  =  [Ma.Ar28  -  avg(Ma.Ar1)]/28  days 

(3) 

AMa.ArMa  ArMax_112  =  [Ma.Ar112  -  avg(Ma.ArMa  ArMax)]/ 

(112  -  Ma. ArMax)  days  (4) 

The  subscript  refers  to  the  age(s)  over  which  the  rate  was 
determined.  Ma.ArMax  corresponds  to  the  age  at  which 
maximum  Ma.Ar  was  obtained  for  a  given  strain.  This  age 
was  used  as  the  initial  time-point  to  calculate  the  rate  of 
marrow  infilling  during  later  growth.  avg(Tt.ArJ  and  avg- 
(Ma.Aq)  refers  to  the  strain  average  Tt.Ar  and  Ma.Ar 
value  at  time-point  i. 
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Table  1.  Body  Weight  and  Morphological  Traits  for  Femoral  Diaphyses  of  Three  Inbred  Mouse  Strains  Between 

Embryonic  Day  18.5  (E18.5)  and  1  Year  of  Age 


Cortical 

Total  area 

Marrow 

Average  cortical 

Strain  (age) 

N 

Weight  (g) 

area  (mm2) 

(mm2) 

area  (mm2) 

Ja  (mm4) 

thickness  (mm) 

A 

E18.5 

5 

0.11  ±  0.02* 

0.16  ±  0.02*'* 

0.05  ±  0.01*'* 

0.004  ±  0.001* 

0.10  ±  0.02 

Id 

15 

1.3  ±0.1 

0.13  ±  0.02* 

0.30  ±  0.05* 

0.17  ±  0.03 

0.010  ±  0.004* 

0.08  ±  0.01 

7d 

12 

2.8  ±  0.5* 

0.22  ±  0.02 

0.50  ±  0.05*'* 

0.28  ±  0.04*'* 

0.028  ±  0.005 

0.11  ±  0.01* 

14d 

8 

6.0  ±  0.4 

0.23  ±  0.05* 

0.59  ±  0.08*'* 

0.35  ±  0.05*-* 

0.035  ±  0.011*-* 

0.10  ±  0.02* 

28d 

8 

10.8  ±  1.6 

0.37  ±  0.03* 

0.82  ±  0.06*'* 

0.45  ±  0.04*'* 

0.076  ±  0.012*'* 

0.13  ±0.01* 

56d 

9 

16.4  ±  2.3 

0.51  ±  0.09* 

0.97  ±  0.10*'* 

0.46  ±  0.03* 

0.119  ±  0.027*'* 

0.18  ±  0.03 

112d 

13 

22.9  ±  3.0+ 

0.74  ±  0.04* 

1.15  ±0.8*-* 

0.41  ±  0.05*-* 

0.189  ±  0.023*'* 

0.25  ±  0.02*'* 

182d 

12 

23.2  ±  1.6™ 

0.79  ±  0.04* 

1.25  ±0.06*'* 

0.46  ±  0.03*'* 

0.224  ±  0.024*'* 

0.26  ±  0.02*'* 

365d 

B6 

E18.5 

9 

24.7  ±  2.4 

0.79  ±  0.06* 

1.28  ±0.11*'* 

0.49  ±  0.07* 

0.228  ±  0.037*-* 

0.25  ±  0.02* 

10 

_ 

0.17  ±  0.04*-* 

0.33  ±  0.04*-* 

0.16  ±  0.04*'* 

0.014  ±  0.004*'* 

0.10  ±  0.03 

Id 

15 

1.3  ±0.1 

0.16  ±  0.03* 

0.35  ±  0.03* 

0.19  ±  0.04 

0.014  ±  0.003* 

0.09  ±  0.02 

7d 

12 

3.1  ±  0.3* 

0.19  ±  0.05 

0.61  ±  0.05* 

0.41  ±  0.05*'* 

0.032  ±  0.008 

0.08  ±  0.02* 

14d 

12 

6.4  ±  0.5* 

0.24  ±  0.04* 

0.90  ±  0.08*-* 

0.66  ±  0.06*'* 

0.061  ±  0.012* 

0.08  ±0.01*-* 

28d 

6 

10.9  ±  1.0 

0.39  ±  0.05 

1.24  ±0.09*'* 

0.85  ±  0.06*'* 

0.133  ±  0.023*'* 

0.11  ±0.01*'* 

56d 

11 

16.8  ±  0.8 

0.57  ±  0.03* 

1.47  ±0.07*-* 

0.90  ±  0.05*'* 

0.226  ±  0.022*'* 

0.16  ±0.01* 

112d 

10 

20.9  ±  0.6* 

0.75  ±  0.03* 

1.55  ±  0.06*'* 

0.80  ±  0.04*'* 

0.297  ±  0.026* 

0.21  ±  0.03*'* 

182d 

9 

21.7  ±  0.9* 

0.76  ±  0.05* 

1.48  ±0.09*-* 

0.72  ±  0.06*'* 

0.283  ±  0.029* 

0.20  ±  0.04*-* 

365d 

10 

25.2  ±  1.8 

0.82  ±  0.06* 

1.82  ±0.10* 

1.00  ±  0.07*'* 

0.387  ±  0.044* 

0.20  ±  0.02* 

C3H 

E18.5 

8 

— 

0.13  ±0.01* 

0.24  ±  0.03*'* 

0.11  ±  0.02*'* 

0.007  ±  0.001* 

0.10  ±0.01 

Id 

15 

1.4  ±0.1 

0.15  ±  0.03 

0.32  ±  0.03 

0.16  ±  0.02 

0.012  ±  0.002 

0.09  ±  0.02 

7d 

9 

4.4  ±  0.5*-* 

0.23  ±  0.06 

0.56  ±  0.07* 

0.33  ±  0.04*'* 

0.032  ±  0.010 

0.10  ±  0.02 

14d 

10 

5.7  ±  0.4* 

0.31  ±  0.03*'* 

0.77  ±0.11*'* 

0.46  ±  0.10*'* 

0.059  ±  0.015* 

0.12  ±0.01* 

28d 

8 

10.4  ±  1.5 

0.43  ±  0.05* 

0.97  ±  0.7*'* 

0.54  ±  0.04*'* 

0.107  ±  0.017*'* 

0.14  ±0.01* 

56d 

10 

15.6  ±2.2 

0.67  ±  0.07*'* 

1.11  ±0.09*'* 

0.44  ±  0.04* 

0.174  ±  0.029*-* 

0.21  ±  0.04* 

112d 

11 

22.4  ±  1.3 

1.08  ±  0.04*  * 

1.41  ±0.06*'* 

0.34  ±  0.03*'* 

0.310  ±  0.025* 

0.34  ±  0.01*'* 

182d 

9 

20.9  ±  0.9* 

1.08  ±  0.05*-* 

1.35  ±  0.05*-* 

0.26  ±  0.04*'* 

0.285  ±  0.024* 

0.37  ±  0.03*-* 

365d 

7 

23.7  ±  1.8 

1.29  ±  0.17*'* 

1.71  ±0.23* 

0.42  ±  0.10* 

0.448  ±0.115* 

0.35  ±  0.04*'* 

*  Significantly  different  from  A,  ANOVA  p  <  0.05. 

*  Significantly  different  from  B6,  ANOVA  p  <  0.05. 

*  Significantly  different  from  C3H,  ANOVA  p  <  0.05. 


Compositional  traits 

Ash  content  was  assessed  for  the  left  femurs  of  all  mice  2 
weeks  of  age  and  older.  These  femurs  had  previously  been 
loaded  to  failure  in  four-point  bending  (data  not  included). 
The  diaphyseal  pieces  were  retrieved  and  cleaned  of  all  soft 
tissues  under  a  stereomicroscope.  The  hydrated,  dried,  and 
ashed  weights  were  determined  as  described  previously/9-* 
Water  content  was  defined  as  the  hydrated  weight  minus 
dried  weight  and  expressed  as  a  percentage  of  hydrated 
weight.  Ash  content  was  determined  as  the  ash  weight  nor¬ 
malized  for  hydrated  weight. 

Statistics 

All  results  are  expressed  as  mean  ±  SD.  Differences  in 
trait  values  among  the  three  genotypes  at  each  time -point 
were  determined  with  a  one-way  ANOVA  and  a  Tukey’s 
posthoc  test.  A  p  value  less  than  0.05  was  considered  sig¬ 
nificant  (GraphPad  Software,  San  Diego,  CA,  USA). 

RESULTS 

Overall  changes  in  femoral  diaphyseal  structure 

All  three  strains  underwent  qualitatively  similar  changes 
in  bone  architecture,  size,  and  shape  of  the  mid-diaphyseal 


cortex  as  a  function  of  age  (Fig.  1).  For  example,  each  strain 
exhibited  an  entirely  woven  or  porous  tissue  microstructure 
at  E18.5.  Flowever,  by  14  days  of  age,  the  woven  tissue  was 
located  only  in  the  posterior-lateral  quadrant,  and  by  28 
days  of  age,  the  entire  mid-diaphysis  exhibited  a  solid,  com¬ 
pact,  lamellar  tissue  structure.  Flowever,  despite  similar 
overall  growth  patterns,  significant  variation  in  bone  size, 
shape,  and  composition  were  observed  among  the  strains 
throughout  growth  (Table  1). 

Prenatal  bone  growth 

Significant  differences  in  total  area  (Tt.Ar),  cortical  area 
(Ct.Ar),  and  polar  moment  of  inertia  (JD)  were  observed 
as  early  as  E18.5  (Fig.  2A;  ANOVA,  p  <  0.05).  B6  mice 
showed  the  largest  Ct.Ar,  Tt.Ar,  and  JD,  whereas  A  mice 
showed  the  smallest.  C3H  mice  more  closely  resembled 
A  mice  at  E18.5.  Between  E18.5  and  1  day  of  age  none  of 
the  morphological  traits  changed  in  B6  mice.  In  contrast, 
A  and  C3FI  mice  showed  increases  in  all  three  traits,  so  that 
by  1  day  of  age  only  A  mice  showed  significantly  smaller 
Tt.Ar,  Ct.Ar,  and  JD  values  compared  with  B6  (Fig.  2B; 
p  <  0.05). 
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equivalent  locations  within  the  femoral  diaphysis  taken  immediately  distal  to  the  third  trochanter.  Note  change  from  a  porous,  woven 
structure  to  a  compact,  lamellar  structure  between  PI  and  P28  days  of  age.  Scale  bar  shown  applies  to  all  images. 
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FIG.  2.  Variation  in  femoral  cortical  area  (Ct.Ar),  total  area 
(Tt.Ar),  and  polar  moment  of  inertia  (JD)  among  A,  B6,  and  C3H 
mice  at  (A)  embryonic  day  18.5  (E18.5)  and  (B)  postnatal  day  1. 
a,  b,  and  c  indicates  a  significant  difference  from  A,  B6,  and  C3H 
mice,  respectively  (ANOVA,  p  <  0.05). 


AGE  (DAYS) 

FIG.  3.  Variation  in  (A)  cortical  area  (Ct.Ar),  (B)  polar  moment 
of  inertia  (JD),  and  (C)  ash  content  across  development  for  A,  B6, 
and  C3H  femurs.  Significant  differences  among  strains  in  Ct.Ar, 
J0,  and  ash  content  are  indicated  in  Table  1. 


Postnatal  changes  in  Ct.Ar  and  J 0 

After  birth,  Ct.Ar  and  JQ  increased  nonlinearly  to  112 
days  of  age  (Figs.  3A  and  3B).  Significant  differences  in 
Ct.Ar  and  J0  values  were  observed  among  the  strains 


throughout  development  (Table  1);  however,  these  differ¬ 
ences  became  more  prominent  at  56  days  of  age.  For  Ct.Ar, 
A  and  B6  mice  showed  similar  values  beyond  7  days  of  age, 
whereas  C3H  mice  showed  larger  Ct.Ar  than  both  A  and 
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FIG.  4.  Variation  in  femoral  total  area 
(Tt.Ar)  and  marrow  area  (Ma.Ar)  among  A, 
B6,  and  C3H  female  mouse  strains.  (A) 
Changes  in  total  area  across  development. 
(B)  Rate  of  change  in  total  area  during  early 
growth  (P1-P28  days  of  age)  and  later 
growth  (28-112  days  of  age).  (C)  Changes  in 
marrow  area  across  development.  (D)  Rate 
of  change  in  marrow  area  during  early 
growth  (P1-P28  days  of  age)  and  growth  af¬ 
ter  maximum  marrow  area  was  achieved 
(Ma.ArMax-112  days  of  age).  Significant  dif¬ 
ferences  among  strains  in  Tt.Ar  and  Ma.Ar 
are  indicated  in  Table  1.  a,  b,  and  c  indicates 
a  significant  difference  from  A.  B6,  and  C3H 
mice,  respectively,  at  the  indicated  time- 
point. 


B6  after  14  days  of  age  (ANOVA,  p  <  0.05).  For  Jo;  B6  and 
C3H  mice  showed  similar  values  beyond  56  days  of  age, 
whereas  A  mice  showed  smaller  values  compared  with  B6 
and  C3H  beyond  14  days  of  age  (p  <  0.05).  These  differ¬ 
ences  in  physical  bone  traits  were  accompanied  by  nearly 
identical  age-related  increases  in  body  weight  among  all 
strains  (Table  1). 

Postnatal  expansion  of  the  periosteal  surface 

Tt.Ar  increased  nonuniformly  with  age  in  a  genotype- 
dependent  manner  (Fig.  4A).  Each  strain  achieved  80-90% 
of  adult  Tt.Ar  and  periosteal  diameter  by  28  days  of  age.  In 
contrast,  all  three  strains  achieved  only  40-50%  of  adult 
weight,  Ct.Ar,  and  JD  by  this  age.  These  data  indicated  that 
the  periosteum  expanded  rapidly  for  all  strains  early  in  life. 
Flowever,  the  rate  of  expansion  at  the  periosteum  (i.e.,  in¬ 
crease  in  Tt.Ar  per  day)  during  early  growth,  and  the  size  of 
Tt.Ar  at  28  days  of  age  differed  significantly  among  the 
strains  (B6  >  C3H  >  A;  ANOVA,  p  <  0.05;  Fig.  4B;  Table 
1).  Between  28  and  112  days  of  age,  the  average  rate  of 
increase  in  Tt.Ar  was  reduced  by  -80%  compared  with  the 
rates  observed  between  PI  and  P28  days  of  age.  During  the 
later  phase  of  growth  (28-112  days),  the  rate  of  Tt.Ar  ex¬ 
pansion  was  greatest  in  C3FI  mice,  whereas  A  and  B6  mice 
showed  similarly  lower  rates  of  expansion  (p  <  0.05).  Age- 
related  changes  in  periosteal  diameter  mirrored  the 
changes  in  Tt.Ar,  as  expected  (data  not  shown). 

Postnatal  expansion  of  the  endosteal  surface 

Like  Tt.Ar,  marrow  area  (Ma.Ar)  also  changed  nonuni¬ 
formly  with  age  and  in  a  genotype-dependent  manner  (Fig. 
4C).  All  mouse  strains  exhibited  similar  Ma.Ar  values  at 
birth  and  exhibited  rapid  postnatal  marrow  expansion  be¬ 
tween  PI  and  P28  days  of  age.  During  this  early  growth 
phase,  B6  femurs  exhibited  the  greatest  rate  of  endosteal 
expansion  (resorption),  followed  by  C3F1  mice,  and  then  A 


mice,  resulting  in  significant  differences  in  Ma.Ar  size  at  28 
days  of  age  (A  <  B6  <  C3H;  ANOVA,  p  <  0.05;  Fig.  4D). 
Maximum  Ma.Ar  values  were  achieved  for  A  and  C3H 
mice  at  -28  days  of  age,  whereas  this  did  not  occur  until  56 
days  of  age  for  B6  mice.  After  reaching  maximal  Ma.Ar 
values,  all  strains  exhibited  decreases  in  Ma.Ar,  indicating 
that  the  endosteal  surface  underwent  a  reversal  in  biologi¬ 
cal  activity  from  net  resorption  to  net  apposition  (i.e.,  in¬ 
filling).  During  the  period  between  the  time  of  maximal 
Ma.Ar  and  112  days  of  age,  the  three  strains  of  mice 
showed  significantly  different  rates  of  endosteal  contraction 
(apposition;  A  <  B6  <  C3FI;  p  <  0.05).  As  a  result,  adult 
Ma.Ar  at  112  days  of  age  was  largest  in  B6  mice  because  of 
the  early  rapid  expansion  of  the  marrow  cavity  and  smallest 
in  C3H  mice  because  of  its  slower  marrow  excavation  and 
larger  fractional  infilling.  Age-related  changes  in  endosteal 
diameter  were  similar  to  those  observed  for  Ma.Ar,  as  ex¬ 
pected  (data  not  shown). 

Age-related  changes  in  cortical  thickness 

Changes  in  cortical  thickness  (Ct.Th)  reflected  the  rela¬ 
tive  rate  of  expansion  of  the  periosteal  and  endosteal  sur¬ 
faces.  The  average  cortical  thickness  increased  nonlinearly 
between  birth  and  1  year  of  age  in  all  three  strains  (Fig.  5). 
Between  E18.5  and  7  days  of  age,  the  average  cortical  thick¬ 
ness  did  not  change  significantly  for  any  of  the  strains  in¬ 
dicating  that  expansion  of  the  marrow  cavity  kept  pace  with 
expansion  of  the  periosteum.  However,  beyond  14  days  of 
age,  cortical  thickness  increased  linearly  until  peak  values 
were  obtained  at  112  days  of  age.  After  14  days  of  age,  C3H 
femurs  exhibited  cortices  with  the  largest  average  thickness, 
A  femurs  exhibited  intermediate  cortical  thickness,  and  B6 
the  smallest  (ANOVA,  p  <  0.05). 

Postnatal  changes  in  ash  content 

In  all  three  strains,  ash  content  increased  with  age  to  112 
days  (Fig.  6).  Genotype-specific  variation  in  ash  content 
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AGE  (DAYS) 

FIG.  5.  Variation  in  femoral  cortical  thickness  (Ct.Th)  among  A, 
B6,  and  C3H  mice  across  development.  Significant  differences  in 
Ct.Th  among  strains  are  indicated  in  Table  1. 
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FIG.  6.  Variation  in  femoral  ash  content  across  development  for 
A,  B6,  and  C3H  femurs.  *B6  <  A  (ANOVA;  p  <  0.05);  **B6  < 
C3H  (ANOVA;  p  <  0.05). 


was  observed  at  the  earliest  time -point  assayed  (2  weeks  of 
age),  and  these  relationships  persisted  throughout  the  re¬ 
mainder  of  the  study. 


DISCUSSION 

This  study  showed  that  genetic  variability  in  inbred 
mouse  adult  bone  traits  (Ct.Ar,  Jo;  and  ash  content)  results 
directly  from  genotype-specific  differences  in  periosteal  and 
endosteal  expansion  and  whole  bone  mineral  accrual,  as 
hypothesized.  Previous  studies  comparing  B6  and  C3H 
mouse  growth  suggested  that  genetic  variation  in  adult 
bone  mass  could  be  explained  by  prepubertal  and  pubertal 
growth  patterns(27)  and  by  different  rates  of  bone  formation 
during  adolescent  growth. <16)  Additional  studies  examined 
B6  skeletal  development  alone. <28,29)  However,  these  stud¬ 
ies  did  not  include  prenatal  and  immediate  postnatal 
growth  (<7  days).  Therefore,  it  remained  unclear  how  in 
utero  growth,  early  bone  cell  (osteoblast  and  osteoclast) 
activities,  and  peripubertal/adolescent  growth  contributed 
to  genetic  variability  in  adult  trait  magnitudes.  By  looking 
across  growth,  including  several  time-points  not  previously 
studied,  we  were  able  to  connect  early  and  later  bone  de¬ 
velopment.  Thus,  this  work  filled  in  gaps  from  previous 


studies  and  provided  new  insight  into  the  genetic  regulation 
of  adult  bone  traits.  Furthermore,  we  studied  both  the  spa¬ 
tial  distribution  of  bone  (JD)  and  matrix  mineralization  (ash 
content).  Previous  work  from  our  laboratory  has  shown 
that  these  traits  in  combination  with  tissue  amount  (Ct.Ar) 
are  important  in  the  determination  of  whole  bone  strength 
and  brittleness,  respectively/9)  Last,  we  included  a  third 
strain  (A),  which  was  previously  shown  to  have  a  small 
adult  bone  morphology  and  high  adult  ash  content/9,14)  to 
further  generalize  how  genetic  background  influences  adult 
bone  traits. 

Our  examination  of  bone  morphology  from  E18.5  to  365 
days  of  age  revealed  that  A,  B6,  and  C3H  inbred  mouse 
strains  showed  qualitatively  similar  transverse  femoral 
growth  patterns.  The  early  rapid  expansion  of  the  periosteal 
and  endosteal  surfaces  followed  by  marrow  infilling  ob¬ 
served  for  the  mouse  was  consistent  with  that  observed  for 
human  long  bones. (21>  As  such,  studying  the  genetic  regu¬ 
latory  pathways  controlling  murine  long  bone  growth  and 
the  attainment  of  peak  bone  properties  should  prove  valu¬ 
able  in  understanding  similar  processes  in  the  human  skel¬ 
eton. 

The  attainment  of  a  particular  combination  of  Ct.Ar  and 
J0  in  adult  inbred  mice  can  be  explained  in  mathematical 
terms  based  on  variability  in  the  relative  expansion  of  the 
periosteal  and  endosteal  surfaces  during  growth  when  the 
femoral  diaphysis  is  modeled  as  a  simple  hollow  cylin¬ 
der/20  22)  Because  Ct.Ar  and  JD  are  related  to  periosteal 
and  endosteal  diameters  differently, 

Ct.Ar  =  Tl.Ar  -  Ma.Ar  oc  DiaPeri2  -  DiaEndo2  (5) 

JD  oc  Tt.Ar2  -  Ma.Ar2  oc  DiaPeri4  -  DiaEndo4  (6) 

cylinders  of  different  diameters  will  necessarily  have  differ¬ 
ent  combinations  of  Ct.Ar  and  JQ  values.  Thus,  variability  in 
adult  Ct.Ar  and  JD  arise  during  growth  necessarily  because 
differences  in  the  relative  expansion  of  the  periosteal  and 
endosteal  surface  (i.e.,  growth  patterns)  lead  to  different 
periosteal  and  endosteal  diameters. 

We  found  that  the  unique  pattern  of  periosteal  and  end¬ 
osteal  growth  exhibited  throughout  life  by  each  genotype 
could  indeed  explain,  in  mathematical  terms,  how  variabil¬ 
ity  in  adult  morphological  trait  combinations  (Ct.Ar  and  JD) 
arises.  For  B6  mice,  the  large  J0  and  small  Ct.Ar  observed 
in  adulthood  resulted  primarily  from  a  large  periosteal  and 
endosteal  expansion  during  early  growth  (before  puberty). 
For  C3H  mice,  an  intermediate  rate  of  periosteal  and  end¬ 
osteal  expansion  during  early  growth  (compared  with  A 
and  B6)  and  a  significantly  greater  periosteal  expansion  and 
marrow  infilling  after  puberty  led  to  the  large  Ct.Ar  and  JD 
values  of  adult  C3H  mice.  As  a  result,  Ct.Ar  and  JD  in¬ 
creased  at  a  faster  rate  in  C3H  mice  and  led  to  a  higher 
Ct.Ar  relative  to  JD  because  a  larger  proportion  of  new 
bone  was  added,  in  a  less  mechanically  efficient  manner,  to 
the  endosteum.  It  is  noteworthy  that  B6  and  C3H  mice 
achieved  similar  adult  JQ  values.  However,  this  occurred 
through  very  different  growth  patterns  and,  presumably, 
different  biological  control  mechanisms.  For  A  mice,  the 
small  adult  Ct.Ar  and  JD  were  a  result  of  having  the  smallest 
periosteal  and  endosteal  expansion  both  before  and  after 
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puberty.  Interestingly,  A  and  B6  mice  achieved  similar 
adult  Ct.Ar  values,  but  this  occurred  through  very  different 
growth  patterns,  suggesting  that  A  and  B6  mice  achieve 
similar  adult  Ct.Ar  through  different  biological  control 
mechanisms.  Together  these  results  show  that  inbred  mice 
build  morphologically  distinct  femurs  by  regulating  bone 
surface  activity  in  a  genotype-specific  manner. 

From  these  data,  we  identified  six  bone  traits  that  helped 
characterize  the  variability  in  transverse  femoral  develop¬ 
ment  among  three  strains  of  inbred  mice  and  that  provided 
additional  insight  into  the  biological  control  mechanisms 
underlying  variability  in  adult  bone  traits.  The  six  traits 
included  (1)  the  initial  size  of  the  bone  at  1  day  of  age,  (2) 
periosteal  expansion  during  the  early  (prepubertal),  rapid 
growth  phase,  (3)  endosteal  expansion  during  the  early 
growth  phase,  (4)  periosteal  expansion  during  the  later 
(peri-  to  postpubertal)  growth  phase,  (5)  marrow  infilling 
during  the  later  growth  phase,  and  (6)  ash  content  (material 
composition)  throughout  growth.  Interstrain  variability  in 
the  first  five  indices  reflects  differences  in  bone  cell  activity 
during  different  growth  phases,  and  the  last  index  reflects 
variation  in  the  tissue  mineralization  process.  Together, 
these  bone  traits  represent  phenotypic  indices  that  will  be 
useful  for  identifying  the  biological  (and  genetic)  mecha¬ 
nisms  that  influence  adult  whole  bone  mechanical  proper¬ 
ties. 

Variability  in  bone  size  and  shape  (Ct.Ar,  JD,  and  TAr)  at 
E18.5  and  1  day  of  age  show  that  these  inbred  mice  begin 
postnatal  skeletal  development  with  different  initial  mor¬ 
phologies.  These  results  suggest  that  bone  growth  and  adult 
bone  mass  might  be  related  to  allelic  variability  in  the  fac¬ 
tors  involved  in  in  utero  development.  Differences  in  tissue 
patterning  and  condensation  have  been  observed  in  the  cra- 
nio-facial  development  of  B6  and  C3H  mice/305  and  similar 
variability  in  tissue  patterning  and  condensation,  cartilage 
formation,  or  initial  bone  formation  could  influence  long 
bone  morphology  at  birth.  Research  in  humans  has  identi¬ 
fied  correlations  between  in  utero  growth,  postnatal 
growth,  and  adult  bone  mass  that  support  the  idea  that 
pre-  and  postnatal  growth  influence  adult  skeletal  proper¬ 
ties. (31_33)  Additional  studies  in  mice  may  be  useful  in  es¬ 
tablishing  the  relationship  between  prenatal  morphology 
and  growth  and  adult  bone  trait  values. 

Investigation  of  early  growth  (1-28  days  of  age)  revealed 
that  prepubertal  growth  played  a  critical  role  in  defining 
adult  trait  values.  Strain  specific  differences  in  adult  exter¬ 
nal  size  were  established  during  early  growth,  suggesting 
that  variation  in  Tt.Ar,  bone  width,  and  maximal  Ma.Ar  are 
likely  the  result  of  allelic  variations  influencing  the  pro¬ 
cesses  controlling  rapid  surface  expansion  before  puberty. 
Measures  of  adult  bone  morphology,  which  have  been  cor¬ 
related  with  fracture  risk/3,345  are  also  largely  determined 
before  puberty  in  humans/215  Therefore,  the  mouse  may 
provide  a  useful  model  for  studying  how  fracture  risk  is 
associated  with  the  biology  of  early  bone  growth  and  the 
establishment  of  external  size. 

Variability  in  ash  content  was  also  established  at  the  ear¬ 
liest  time-point  studied  (2  weeks  of  age)  and  was  consistent 
with  the  variability  observed  in  adulthood.  This  suggested 
that  differences  in  mineralization  and  material  composition 


among  the  strains  might  be  observed  even  earlier  in  life, 
possibly  before  birth,  or  within  the  first  few  days  of  post¬ 
natal  growth.  We  previously  showed  that  a  4.0%  variation 
in  ash  content  among  inbred  strains  was  correlated  with  a 
-3.4-  and  -1. 9-fold  variation  in  whole  bone  ductility  and 
toughness,  respectively. (14>  Therefore,  the  consistently 
lower  B6  femoral  ash  content,  compared  with  A  and  C3H 
mice,  implied  that  B6  femurs  were  constructed  of  bone  tis¬ 
sue  with  a  lower  material  stiffness  but  greater  ductility5355 
throughout  life.  Together,  these  results  suggest  that  differ¬ 
ent  strains  of  inbred  mice  might  harbor  allelic  variations  in 
the  genetic  factors  inherent  to  or  governing  mineral  accrual. 

The  rate  and  extent  of  periosteal  and  endosteal  expan¬ 
sion  during  early  growth  appeared  to  be  well  matched  to 
ash  content  values,  suggesting  that  material  composition 
may  be  functionally  coupled  to  bone  growth  patterns.  B6 
mice,  with  the  most  rapid  prepubertal  surface  expansions 
and  largest  Tt.Ar  and  J0  values  at  28  days  of  age,  showed 
small  ash  values.  In  contrast,  the  slower  expanding  bones  of 
A  and  C3H  mice  showed  larger  ash  values.  From  a  me¬ 
chanical  basis,  these  results  are  not  entirely  surprising  be¬ 
cause  the  addition  of  an  equivalent  amount  of  tissue  on  the 
outer  surface  of  a  large  bone  will  result  in  a  greater  increase 
in  JD,  and  therefore  bone  stiffness/29,355  compared  with  a 
similar  addition  on  a  smaller  bone.  Thus,  the  rapid  expan¬ 
sion  of  the  periosteum  and  endosteum  in  B6  mice  estab¬ 
lished  a  more  efficient  growth  pattern  compared  with  A 
and  C3H  mice.  This  greater  structural  efficiency  would  sub¬ 
stantially  increase  whole  bone  stiffness  and  presumably  af¬ 
ford  B6  mice  the  ability  to  reach  a  mechanically  functional 
end  state  with  material  of  reduced  ash  content  and  thus 
reduced  tissue  level  stiffness. 

The  variability  in  femoral  growth  patterns  observed  after 
the  onset  of  puberty  were  entirely  consistent  with  previ¬ 
ously  reported  data/16,27,295  Furthermore,  our  comparison 
of  early  and  later  growth  suggests  that  female  C3H  mice 
might  harbor  allelic  variations  that  significantly  alter  their 
growth  patterns  after  puberty  and  predispose  them  to  dras¬ 
tically  increased  bone  growth  in  association  with  sexual 
maturation.  Additional  studies  are  required  to  test  if  pu¬ 
bertal  maturation  and  sex  hormones  differentially  affect  the 
biological  processes  of  skeletal  growth  and  contribute  to  the 
variability  in  adult  bone  traits  among  these  strains. 

Significant  variation  in  the  aging  of  the  inbred  mouse 
skeleton  was  also  observed.  A  mice  showed  no  changes  in 
Tt.Ar  and  Ma.Ar  values  between  182  and  365  days  of  age, 
suggesting  little  periosteal  expansion  and  endosteal  expan¬ 
sion  during  this  time,  whereas  B6  and  C3H  mice  showed 
large  increases  in  Tt.Ar  and  Ma.Ar.  The  age-related  expan¬ 
sion  of  B6  and  C3H  femurs  is  similar  to  that  observed  for 
the  human  skeleton/21,365  These  results  suggest  that  the 
biological  processes  associated  with  skeletal  aging  vary  with 
genetic  background. 

Although  female  A,  B6,  and  C3H  mice  exhibited  a  wide 
range  of  skeletal  phenotypes,  further  studies  are  warranted 
to  test  how  growth  patterns  vary  among  other  strains  of 
mice,  how  growth  patterns  vary  between  male  and  fe¬ 
males,' (21'23’24)  and  how  each  phase  of  growth  is  related  to 
biomechanical  properties.  In  a  standard  cylindrical  model 
of  growth,  measures  of  Tt.Ar  and  Ma.Ar  are  sufficient  to 
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describe  how  a  single  species  of  bone  cells  (osteoblasts  or 
osteoclast)  influence  surface  changes.  However,  the  trans¬ 
verse  growth  of  the  femoral  mid-diaphysis  occurs  through 
both  surface  expansion  and  cortical  drift/37’38-1  In  regions  of 
long  bones  exhibiting  drift,  both  apposition  and  resorption 
occur  at  different  locations  along  each  surface;  therefore, 
the  relative  contribution  of  osteoblasts  and  osteoclasts 
working  on  both  the  periosteum  and  endosteum  must  be 
considered  in  the  analysis  of  Tt. Ar  and  Ma.Ar,  respectively. 

From  this  study,  we  conclude  that  different  strains  of 
inbred  mice  build  distinct,  yet  mechanically  functional,  fe¬ 
murs  by  regulating  bone  surface  activity  and  mineral  ac¬ 
crual  in  a  genotype-specific  manner.  These  genotype- 
specific  growth  patterns  led  directly  to  particular 
combinations  of  adult  bone  traits,  and  thus,  specific  reper¬ 
toires  of  whole  bone  mechanical  properties.  By  looking 
across  all  ages  of  growth,  from  prenatal  development  to 
adulthood,  we  identified  six  bone  traits  that  can  be  used  to 
explain,  in  a  systematic  manner,  the  differences  in  femoral 
growth  patterns  among  inbred  mice.  This  characterization 
of  femoral  growth  for  A,  B6,  and  C3H  mice  throughout  life 
revealed  how  and  when  specific  combinations  of  adult  bone 
traits  arise  in  genetically  distinct  animals.  Thus,  the  inbred 
mouse  provides  a  valuable  model  to  determine  how  genetic 
background  influences  the  biological  processes  involved  in 
establishing  variability  in  adult  bone  traits,  and  thus,  whole 
bone  mechanical  properties. 
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